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Assembly  and  Application  of  a  Shock  Data  Analysis  System 

1.  Introduction 

Measured  data  from  shipboard  shock  environments  usually  consist  of 
analogue  magnetic-tape  recordings  of  acceleration  and  strain  versus  time. 
Usually,  little  analysis  is  required  of  strain  recordings  save  spectral  analy¬ 
sis  and  perhaps  sums  and  differences.  From  the  accelerometer  recordings, 
however,  it  is  usually  necessary  to  extract  velocities  and  displacements,  and 
often  shock  spectra,  as  well.  Shock  data  analysis  facilities  tend  to  run  to 
two  extremes  of  size  and  capacity.  At  one  end  is  the  large  digital  processing 
facility  where  hundreds  of  records  (from  an  at-sea  ship  shock  test,  for 
example)  are  processed  in  batch,  usually  after  the  test  is  completed.  At  the 
under  end  is  the  small,  often  analogue,  facility  which  can  do  limited  analysis 
of  a  small  number  of  records,  although  perhaps  very  slowly.  The  latter  may  be 
used  to  advantage  for  analyzing  a  few  selected  critical  records  in  the  field, 
or  in  the  shock  simulation  lab. 

A  need  exists  for  a  capability  in  the  middle  ground  -  a  system  capable  of 
doing  analysis  on  a  modest  number  of  records  in  a  reasonable  time,  which 
accepts  analogue  magnetic  tapes  as  input  and  produces  report-ready  graphs  and 
tables  as  output,  and  which  is  field-transportable  if  necessary. 

This  report  describes  such  a  system  which  has  been  assembled  with 
standard,  off-the-shelf  instruments  which  are  marketed  nationwide  and  easily 
available.  It  is  not  in  any  way  optimized  for  the  specific  task  of  shock  data 
analysis,  hence  may  not  be  competitive  in  speed,  cost,  etc.  to  a  system  which 
has  been  optimized  in  hardware  and/or  software  to  this  specific  application. 

It  furnishes  an  example  of  what  can  be  done  along  these  lines  with  virtually 
no  expenditure  of  R  &  D  effort,  and  is  an  example  in  that  any  of  its  component 
instruments  could  be  replaced  by  a  different  supplier's  product  with  little 
change  in  operation  or  performance. 

2.  Background 

2.1.  Shock  motions  may  be  described  in  terms  of  acceleration  -, 
velocity  -,  or  displacement-time  histories,  and  in  all  three  parameters  may 
vary  greatly  in  waveform  complexity.  Simpler  shock  motions  may  be  closely 
defined  by  a  few  waveform  parameters,  and  these  may  also  provide  a  reliable 
basis  for  comparison  of  shock  motions.  Shipboard  shock,  however,  is  extremely 
complex.  Motions  are  typified  by  waveforms  of  high  frequency  range,  high 
dynamic  range  and  long  duration  which  cannot  be  defined  adequately  by  a  few 
parameters.  The  additional  descriptions  found  most  useful  for  a  motion  of 
this  type  are  the  shock  spectrum,  which  may  be  interpreted  as  a  measure  of  the 
effects  of  a  shock  motion,  and  the  design  shock  spectrum,  which  provides  the 
data  needed  to  design  a  structure  to  withstand  it.  1!ach  provides  a  basis  for 
comparing  shock  severities,  and  applies  equally  well  to  shock  motions  of  all 
degrees  of  complexity  of  waveform.  The  general  nature  of  shock  motions,  shock 
spectra  and  design  shock  spectra  are  discussed  in  Refs.  (1)  and  (2). 

2.2.  While  a  shock  motion  may  be  identified  equally  well  by  its  acceler- 
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ation,  velocity  or  displacement  waveform,  one  of  them  will  usually  be  most 
convenient  for  a  given  application.  However,  in  most  situations  only  acceler¬ 
ation  can  be  measured  properly.  This  is  almost  always  the  situation  in  the 
case  of  shipboard  shock,  and  the  bulk  of  data  in  this  field  will  consist  of 
records  of  acceleration  and  strain  versus  time,  with  perhaps  a  small  fraction 
being  records  of  velocity  versus  time.  Thus,  an  analysis  system  for  shipboard 
shock  motion  records  must  be  able  to  produce  time  histories  of  acceleration, 
velocity  and  displacement,  as  well  as  shock  spectra  and  design  shock  spectra, 
given  input  signals  representing  acceleration  -  or  velocity-time  histories, 
and  preferably  from  displacement-time  histories  as  well.  Since  shipboard 
shock  motions  are  among  the  most  demanding  in  terms  of  functional 
requirements,  an  analysis  system  capable  of  handling  them  will  also  be  able  to 
handle  shock  motions  of  most  other  tjrpes. 

3.  Analysis  System  Characteristics 

3.1.  The  principal  application  of  the  system  is  the  analysis  of  shock 
motions  measured  on  laboratory  shock  machines.  This  includes  those  which 
generate  simple-pulse  shock,  but  most  prominent  are  those  which  generate 
shipboard  shock  type  environments:  the  Lightweight  Shock  Machine  (LWSM),  the 
Medlumwelght  Shock  Machine  (MWSM) ,  the  Floating  Shock  Platform  (FSP)  and  the 
Large  Floating  Shock  Platform  (LFSP)  (Ref.  3).  Tests  on  these  machines  will 
be  lightly  instrumented  compared  to  at-sea  tests  of  ships  and  submarines,  or 
some  of  the  more  specialized  test  vehicles.  The  typical  instrumentation  suite 
may  be  a  few  channels  on  the  slmpler-pulse  machines  and  the  LWSM,  up  to  a 
dozen  on  the  MWSM  and  FSP  and  a  few  dozen  on  the  LFSP,  while  on  some  occasions 
several  times  as  many  channels  may  be  required.  Normally,  about  one-third  of 
the  total  channels  will  fall  in  the  category  of  equipment  input  motions,  while 
most  of  the  remainder  will  be  equipment  response  motions  and  structural 
strains. 


3.2.  The  function  of  the  analysis  system  is  to  accept  analogue  input 
signals,  perform  whatever  processing  is  required,  and  produce  output  represen¬ 
tations  of  the  desired  quantities,  scaled  in  the  appropriate  engineering 
units,  in  graphical  or  tabular  form,  and  in  video  or  hard-copy  display.  All 
signals  require  A-D  conversion  and  storage,  scaling  and  display.  All  motion 
signals  require  in  addition  computation  and  display  of  acceleration,  velocity 
and  displacement,  each  properly  scaled,  which  requires  the  functions  of 
integration  and  possibly  differentiation.  Input  motion  signals,  and  some 
response  motion  signals,  also  require  the  computation  of  maximax  and  residual 
shock  spectra.  This  last  computation  is  by  far  the  most  complex  and  time- 
consuming,  but  is  needed  for  only  a  few  channels.  For  a  suite  of  a  dozen 
channels  of  assorted  information,  the  system  will  perform  a  complete  analysis 
in  an  hour  or  two. 


Shock  spectrum  analyzers  are  available  commercially,  and  are  typically 
compact  units  and  very  fast.  As  stand-alone  systems,  they  are  lacking  in  that 
they  cannot  do  a  complete  analysis  (integration  and  scaling,  etc.)  and  have  no 
display.  As  components  for  calculating  shock  spectra  they  are  wanting  in  that 
they  cannot  compute  undamped  spectra  and  their  analysis  frequencies  are  fixed. 
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3.3.  Since  Che  analysis  system  is  an  assembly  of  component  units  designed 
for  use  as  independent,  general  purpose  instruments,  it  has  a  great  degree  of 
flexibility.  With  different  switch  settings,  it  would  function  equally  well 
for  analyzing  earthquake  motions  or  blast  loadings.  A  brief  description  of 
each  functional  block  of  the  system  and  its  capabilities  follows. 

3.3.1.  Magnetic-Tape  Recorder 

The  primitive  form  of  shock  data  records  is  usually  FM 
magnetic  tape.  The  recorder  used  with  the  analysis  system  has  14- 
track  record/reproduce  capability,  plus  edge  tracks  for  voice  annota¬ 
tion.  The  normal  configuration  has  two  housekeeping  channels,  one  of 
them  DR  for  a  tape-speed  servo  reference  signal  and  Che  other  an  FM 
channel  devoted  to  time-code  in  IRIR  B  format,  and  12  FM  data 
channels.  All  FM  channels  are  TRIG  Intermediate  Band  I.  Normal 
recording  speed  is  30  ips,  giving  a  nominal  bandwidth  of  0  -  10 
KHz.  The  recorder  can  be  operated  at  tapes  speeds  from  15/16  to  120 
ips,  allowing  Che  traditional  analysis  tricks  of  data  speed-up  and 
slow-down;  record  density  can  be  changed  by  rearranging  plug-in 
components  on  the  electronics  cards.  Output  level  is  normally  set 
for  1_^  rms  for  ±40%  peak-to-peak  (sinusoidal)  modulation.  This 
particular  unit  was  selected  because  it  performs  a  set  of  canned 
calibration  and  set-up  routines  under  control  of  a  built-in  micropro¬ 
cessor,  and  because  it  is  a  portable  unit.  Both  features  are 
convenient  for  field  operations. 

3.3.2.  Auxiliary  Equipment 

The  principal  auxiliary  equipment  units  are  a  time-code  gener- 
aCor/reader,  a  direct-writing  oscillograph,  and  a  storage  CRO.  In 
generator  mode,  the  time-code  unit  produces  an  IRIG  B  code  for 
recording  in  the  magnetic  tape.  In  reader  mode,  it  monitors  the 
recorded  time-code  and  produces  a  pulse  when  a  chosen  time  value 
occurs.  This  time  can  be  selected  to  the  millisecond  and  provides  a 
start  signal  to  the  transient-storage  unit  (discussed  below)  and 
facilitates  synchronization  of  the  recorded  data  signals.  The  other 
units,  the  CRO  and  oscillograph,  provide  quick-look  display  of  the 
recorded  data  signals  and  time-code.  The  CRO  is  a  two-channel, 
single  time-base  type  similar  to  many  currently  available.  The 
oscillograph  produces  strip  charts  of  up  to  12  channels  with  a  band¬ 
width  of  0  -  5  KHz  at  paper  speeds  from  0.1  to  160  ips. 

3.3.3.  Data  Filter 

The  data  signals  output  from  the  tape  recorder  are  passed 
through  an  8-channel  data  filter.  Each  channel  has  a  low-pass 
characteristic  with  cutoff  frequency  selectable  from  1  Hz  to  15  KHz, 
and  is  direct-coupled.  The  unique  features  of  this  unit  are  that  the 
filter  characteristic  is  optimized  for  linear  phase-shift  within 
pass-band,  minimizing  envelope  distortion,  and  that  phase-shift 
between  channels  is  closely  controlled  at  no  more  than  2®  between  any 
pair.  The  filter  characteristic  also  features  a  cutoff  rate  of  43 
db/octave,  so  that  aliasing  is  negligible  for  sampling  frequencies  of 
4.2  times  the  cutoff  frequency  or  higher.  Each  channel  has  unity 
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gain,  and  can  accept  inputs  of  up  to  10  ^  .  As  used  In  the  analysis 
system,  the  cutoff  frequency  Is  normally  set  to  IKHz. 

3.3.4.  Transient  Storage  Unit 

The  transient-storage  unit  contains  seven  Independent  AD  and 
digital  storage  channels  with  a  common  tine-base.  Each  channel  can 
be  operated  in  an  Inhibit  mode  which  preserves  the  existing  stored 
signal.  The  unit  can  be  operated  in  a  self-triggered  mode,  where  the 
signal  in  any  channel  can  be  used  to  initiate  the  capture  of  the 
signals  in  all  channels  which  are  not  inhibited.  As  used  in  the 
analysis  system,  the  capture-initiate  signal  is  derived  from  the 
time-code  reader  to  preserve  commonality  of  time-base  between  succes¬ 
sive  batches  of  signals.  Upon  initiation,  each  channel  captures  4096 
10-bit  samples  of  its  signal,  for  a  dynamic  range  of  60  db.  As  used 
in  the  analysis  system,  the  channel  full-scale  range  is  ±  2  ^  ;  it 
can  be  varied  from  ±0. 1  to  The  sampling  period  used  is  100 

mlcrosec.,  for  a  total  record  length  of  400  millisec.;  the  sampling 
period  can  be  varied  from  5  mlcrosec.  to  2.5  sec.  Each  channel 
provides  outputs  of  the  digitized  "staircase"  and  smoothed  analogue 
versions  of  the  stored  signal  waveform.  As  used  in  the  analysis 
system,  the  signal  stored  in  any  selected  channel  is  read  out  on  the 
IEEE-488  bus  in  bit-parallel,  character- serial  ASCII  code.  This 
particular  unit  was  selected  because  of  its  7-channel  capacity  and 
the  fact  that  each  A-D  converter  and  storage  channel  is  functionally 
independent  from  the  others. 

3.3.5.  Processor  and  Display  Unit 

The  processor  and  display  unit  is  a  large-screen  desk-top 
computer  and  peripherals.  The  computer  has  64K-byte  memory  and  a 
resident  BASIC  interpreter,  leaving  some  53K-bytes  of  usable 
memory.  Additional  ROM  packs  provide  signal  processing  routines, 
including  Integration  and  differentiation,  and  text-editing 
capability.  The  principal  display  is  a  built-in  large  screen  CRT. 

The  peripherals  include  a  hard-copy  unit  which  gives  a  permanent 
record  of  the  CRT  display  with  reduced  size,  and  a  digital  X-Y 
plotter  which  gives  high  quality  graphs  up  to  11  x  17  in.  Other 
peripherals  are  a  high  capacity  disc- file  manager,  which  allows  a 
large  library  of  programs  to  be  resident  on  floppy  disc,  and  an 
additional  tape-cassette  drive-unit,  allowing  more  flexibility  in 
data  transfer.  All  of  these  components  are  produced  by  the  same 
manufacturer,  and  are  configured  and  interfaced  to  form  a  compatible 
system.  This  (processor  and  display)  unit  was  selected  because  it 
possesses  the  requisite  computational  power  and  has  an  excellent 
repertoire  of  graphics  routines.  Other  features  are  good  man-machine 
interface,  BASIC  interpreter  and  large-screen  display.  These  seem  to 
assist  the  Intended  user  profile,  engineering/technical  personnel 
without  extensive  familarlty  with  computing  systems,  in  learning  to 
operate  this  analysis  system  with  a  minimun  of  Instruction.  The 
existence  of  a  comprehensive  assortment  of  proven  peripherals  was 
also  a  factor  of  consideration.  Another  feature  of  the  computer  is 
that  it  can  also  be  operated  as  a  high-speed  (9600  baud)  intelligent 
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terminal,  should  operation  with  a  larger  host  computer  be 
desirable.  The  last  peripheral  in  the  system  is  a  high-speed 
printing  terminal  from  another  manufacturer.  When  used  with  the 
anal3wlfr .system,  it  functions  as  a  moderately  fast  line-printer. 

3.4.  The  functional  block  diagram  of  the  analysis  system  is  in  shown  in 
Fig.  1. 

4.  Computation  of  Shock  Spectra 

4.1.  Recalling  that  the  shock  spectrum  of  a  shock  motion  may  be 
interpreted  as  the  maximtm  relative  displacements  of  a  set  of  massless  single- 
degree-of-freedom  systems  whose  bases  are  forced  to  follow  the  motion,  it  is 
apparent  that  shock  spectra  may  be  derived  for  any  degree  of  damping.  For  the 
purpose  of  comparing  shock  motion  environments  and  specifying  design  shock 
spectra,  the  undamped  version  is  preferred,  since  the  damped  spectra  are  less 
sensitive  to  variations  in  the  input  motion.  Shock  spectra  can  be  computed 
from  the  differential  equation  of  motion,  or  from  the  corresponding  integral 
equation.  Either  form  can  be  expressed  to  Involve  the  input  motion  in  terns 
of  acceleration,  displacement  or  velocity.  The  procedure  used  here  is  that  of 
numerical  solution  of  the  differential  equation  relating  relative  displacement 
to  input  accleratlon.  The  computational  procedure  requires  the  input  to  the 
program  as  the  acceleration-time  history  of  the  shock  motion.  The 
computational  procedure  uses  the  algorithm  described  in  Ref.  (4), 

The  equation  of  motion  for  an  undamped  oscillator  of  frequency  o  is 
( 1)  X  +  ai*X  =  -Z 


with  the  solution  for  the  relative  displacement  X 
(2)  X(t)  =  Xg  cos  (Dt  + -^  sin  wt  -  —  sin  u)(t  -  t)  dr. 

The  numerical  integration  equation  derived  by  G.  O'Hara  (Ref.  4)  is 


(3) 

=  X^u)  cos  uh  +  Xj^  sin  uh 


C 

n(  1  -  cos  luh)  ..i  /I  +  cos  uh  sin  (oh\ 

Vn - ZSIH - 


with  a  parabolic  approximation  for  the  digitized  input  function 

-  Z„ 


Z  (t); 


and 

(note  that 
so  that 


Sn  -  Sn-1  -  Zn-t-1  "  2Zn  +  K-1 

is  the  second  difference,  not  a  square) 


(4) 


Z(t) 
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In  these  equations,  h  is  the  time  interval  of  integration. 

By  eliminating  in  (3),  R.  Bort's  implementation  of  equation  3  (Ref.  5) 
is 

(5)  X(t  +  h)  =  2  X(t)  cos  (Dh  -  X(t-h)  +  h  [Z  (t  +  h/2)  -  Z(t  -  h/2)] 

Equation  (5)  has  been  applied  in  this  program  because  it  consists  of  only 
three  additions.  As  implemented  here,  the  term  ^n  brackets 
[Z  (t  +  h/2)  -  Z  (t  -  h/2)l  is  identified  as  hZ  (t) ,  a  sample  point  from  the 
digitized  input  acceleration  record.  Before  processing  input  data  the 
2  cos  u)h  terra  is  determined  and  stored  in  array  C.  All  scaling  factors  are 
deferred  to  the  end  of  the  calculation  for  each  analysis  frequency. 

In  order  to  cut  down  calculation  time,  the  fact  that  lower  analysis 
frequencies  require  fewer  periods  to  reach  peak  response  than  higher 
frequencies  is  employed.  The  number  of  periods  per  analysis  frequency 
considered  in  the  processing  routine  is  given  by 

<'■>  «  -  -nlg-F 

where  F  =  (»)/2it  . 

The  second  procedure  to  speed  up  the  program  selects  data  points  with  respect 
to  cutoff  frequency  and  analysis  frequency.  The  calculation  time  has  been 
reduced  this  way  from  about  50  min.  to  8  min.  while  producing  only  a  minor 
error  (see  Fig.  2A,B). 

5.  System  Operator's  Manual 

5.1.1.  It  is  assumed  that  a  time  record  of  the  shock  motion  - 
including  a  calibration  step  -  is  available  on  an  analog  tape  and  a 
time  code  is  recorded  on  a  separate  track  of  the  tape.  Playing  back 
the  record  (or  calibration  step)  and  the  time  code  on  the  Visicorder, 
determine  the  millisecond  preceding  the  starting  time  of  the 
record.  Set  the  time  code  generator  to  open  the  time  gate  and 
trigger  the  transient  recorder.  The  transient  recorder  accomplishes 
the  A/D  -  conversion  and  stores  the  digitized  time  record  in  its 
memory.  An  interface  bus  allows  the  computer  to  have  access  to  the 
seven  channels  (memories)  of  the  transient  recorder.  Further  data 
processing  is  done  under  the  control  of  the  program  listed  in  the 
Appendix. 

5.1.2.  This  program  is  written  in  BASIC  and  runs  on  the  desktop 
computer  Tektronix  4054.  The  three  parts  of  the  program  are  stored 
on  the  disc  in  three  files  (a  back  up  tape  cartridge  contains  the 
program  in  binary  code).  Part  I  reads  the  calibration  step  and  the 
time  record  from  the  memory  of  the  transient  recorder,  and  calculates 
shock  spectra.  Part  II  contains  all  output  procedures.  Part  III  is 
basically  the  same  as  Part  I  except  that  it  reads  a  calibrated 
digital  time  record  from  the  external  tape  unit  rather  than  from  the 
transient  recorder.  Both  Part  I  and  Part  III  use  the  output 
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•  CODE  9*42 
29-r)*lt-SI  12.49 


(a)  Maximax  shock  spectrum 


CODE  SM2 


(b)  Residual  shock  spectrum 

introduced  by  data  decimation.  Maximax  (a)  and  residual  (b)  shock  spectrum  computed  from  the 
complete  (4096)  and  decimated  data  sets  for  a  half-sine  acceleration  pulse. 


procedures  in  Part  IT. 

5.2  Starting  Up 

The  computer  and  at  least  half  of  the  units  connected  vd.th  the 
computer  by  interface  bus  must  be  turned  on.  To  run  the  program,  set  the 
clock  on  the  disc  drive,  mount  the  unit,  load  the  program  and  enter  "RUN". 
Example:  Enter:  CALL  "SETTIM",  "20  -  MAR  -81  □  8:30" 

Press :  RETURN 

Enter;  CALL  "MOUNT",*)  ,  A$ 

Press;  RETURN 

Enter;  OLD  "  0  SPRGM/PART  I" 

Press:  RETURN 
Enter:  RUN 
Press:  RETURN 

Message  on  the  Screen: 

SPRGM/PART  I  is  loaded  (default  conditions)! 

The  corresponding  procedure  for  using  binary  tape  is  either; 

Enter;  FIND  2 
Press;  RETURN 
Enter;  CALL  "BOLD" 

Press:  RETURN 
Enter:  RUN 
Press:  RETURN 
or  simply  press:  AUTOLOAD 

While  running  the  program,  the  computer  indicates  a  need  of  communication 
with  the  user  by  a  "beep"  tone  and  a  message  on  the  screen. 

The  symbol  O  indicates  one  space. 

5.3  Default  Conditions 

Default  conditions  are  listed  in  line  290  to  500  in  the  Appendix.  All 
default  parameters  may  be  changed  at  any  time  by  entering  the  desired  values 
from  the  keyboard  (e.g.  entering  F2  »  500  changes  the  upper  analysis  frequency 
to  500  Hz).  The  current  value  of  every  parameter  can  be  checked  by  entering 
the  parameter  and  pressing  RETURN:  the  current  value  of  the  parameter  appears 
on  the  screen. 

The  most  important  parameters  defined  by  default  conditions  are  explained 
below : 

FI,  F2 :  Lower/upper  analysis  frequency:  can  be  any 

positive  value  between  1  and  1000. 

C0;  Number  of  analysis  frequencies  between  FI  and 

F2:  depends  on  memory  available  (Part  I  up  to 
about  I 10,  Part  III  about  160).  C®  *  1  is 
possible  also. 

F8 ;  Cut  off  frequency  of  low  pass  filter  (default 
1000  Hz). 

Nl,  N2 :  Lower/upper  value  of  jr-axis  in  log  display 
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N6 ,  N7 :  Lower/upper  value  of  x-axis  in  log  display 

D:  Primary  address  for  output  devices: 

32  screen 

3  ext.  tape 

4  plotter 
40  printer 

P:  selects  grid  in  log  display: 

0<P<in  sets  grid  line  at  P 
P>10  no  grid 
(P  must  not  be  zero!) 

0:  indicates  the  part  of  the  program  currently 

loaded  in  the  memory 

1.  Part  I 

2.  Part  II 

3.  Part  III 

(Do  not  change  parameter  0  from  the  keyboard) . 
M9:  Calibration  matrix. 

A2 :  Offset-correction  matrix. 

5.4.  User  Defined  Keys 


The  first  row  of  user  defined  keys  (No.  11  through  15)  are  used  with 
program  Part  I  or  III.  All  other  user  defined  keys  (No.  1  through  10  and  16 
through  20)  belong  to  program  Part  II.  User  Defined  Key  No,  4  and  _5_  allow 
switching  from  program  Part  I  or  III  to  Part  II  or  vice  versa.  Pressing  a 
user  defined  key  which  belongs  to  a  program  part  not  currently  loaded  in  the 
memory  causes  an  error  message  to  appear  on  the  screen. 

The  following  describes  each  user  defined  key  and  the  system's 

reaction. 


5.4.1.  First  Row  (No.  11  through  15) 

User  Defined  Key  No  11  triggers  the  calibration  routine  in  program 
Part  I.  After  the  time-record  has  been  defined  as  a  displacement, 
velocity  or  acceleration  record  there  are  two  ways  to  fill  the 
calibration  matrix  M9  with  the  calibration  factors  for  each  channel. 

First  the  calibration  factor  may  be  entered  from  the  keyboard  in  input 
units  per  1  Volt,  or  secondly,  the  calibration  step  may  be  read  from  each 
channel  of  the  transient  recorder.  In  the  latter  case  the  calibration 
level  has  to  be  entered  in  appropriate  units  before  the  computer  reads 
the  content  of  each  channel  and  displays  it  on  the  screen.  The  user  may 
define  the  calibration  level  by  means  of  setting  the  crosshair  first  to 
the  upper  and  then  to  the  lower  calibration  level.  The  calibration 
factor  (i.e.  input  unit/digit)  for  this  particular  channel  will  be 
assigned  to  an  element  of  calibration  matrix  M9.  The  program  starts  this 
procedure  all  over  again  until  the  calibration  steps  of  all  7  channels 
have  been  read.  If  the  program  encounters  a  channel  set  to  ground  or 
calibration  mode  an  appropriate  message  appears  on  the  screen  and  no 
calibration  factor  will  be  defined  for  this  particular  channel.  The 
program  proceeds  with  the  next  channel. 

User  Defined  Key  No.  12  reads  the  time  record  from  a  defined 
channel  of  the  transient  recorder  and  displays  a  plot  of  acceleration. 
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velocity  and  dlsplacement^vs  time  on  the  screen.  Also  each  maximum  and 
minimum  value  Is  printed.  When  running  program  Part  I,  the  calibration 
routine  (key  11)  must  be  completed  before  using  Key  tto.  12.  With  program 
part  III,  User  Defined  Key  No.  12  causes  the  computer  to  read  the  time 
record  (Including  other  Important  Information  like  calibration  factor, 
time  step,  etc.)  from  the  external  tape. 

The  read  routine  has  to  be  completed  for  each  record  before 
calculating  shock  data  by  pushing  User  Defined  Key  No.  13.  The  computer 
then  calculates  shock  data  of  the  time  record  that  has  been  Input  by 
means  of  Key  12. 

Using  the  mathematical  algorithm  described  above,  the  program 
computes : 

1 .  Analysis  Frequency 

2.  Pseudo  velocity  Maxlmax-Shock  Spectrum 

3.  Rel.  Velocity  -  Residual- Shock  Spectrum  (or 
Fourier  Transform) 

4.  Rel.  Displacements 

5.  Abs.  Acceleration 

6.  Time  when  maxtmm  occurs 

7.  Phase 

While  calculating,  the  computer  displays  on  the  screen  the  serial 
number  of  analysis  frequency  (from  1  at  the  lowest  to  Zt>  at  the 
highest),  the  analysis  frequency,  the  pseudo-velocity,  and  the 
number  of  data  points  used  for  the  calculation  of  shock  spectra  for 
each  analysis  frequency.  The  computing  time  is  also  given  at  the 
end  of  the  procedure.  ^ 

The  number  of  analysis  frequencies  (1  to  160  ),  as  well  as  any 
Interesting  range  of  analysis  frequencies  between  1  and  1000  Hz,  can 
be  chosen.  (Note  that  F2  has  to  be  larger  than  F^). 

The  time  record  and  all  Important  information,  such  as 
calibration  factor,  time  step,  etc.  can  be  stored  on  a  defined  file 
of  the  external  tape  by  pushing  User  Defined  Key  No.  14.  Additional 
analysis  can  be  done  on  this  time  record  by  means  of  part  III  of  the 
program  (refer  to  User  Defined  Key  No.  12).  An  offset  correction  of 
the  time  record  can  be  made  with  the  procedure  of  User  Defined  Key 
No.  15.  The  first  integral  of  the  time  record  is  displayed  on  the 
screen  and  the  user  may  define  the  zero  line  by  means  of  the  cross 
hair.  After  pushing  any  key  the  computer  calculates  an  offset 
correction  for  this  particular  channel  and  considers  that  at  all 
later  operations.  This  procedure  may  be  repeated  several  times. 

*  The  title  or  desc’^lption  of  the  time  record  is  stored  in  string  W$  and 
appears  on  the  top  of  the  display.  It  may  be  changed  anytime  by  redefining  W  $ 

*Depend3  on  memory  available  for  the  calculation  (program  part  III  allows  up 
to  160  analysis  frequencies). 
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However,  the  user  has  to  realize  that  this  procedure  requires  some  actual 
value  to  be  known. 

5.4.2.  Second  Row  (No.  1  through  5) 

The  results  of  the  shock  spectra  computation  (Key  No.  13)  can 
be  saved  by  pushing  User  Defined  Key  No.  1.  This  causes  the 
computer  to  write  all  data  on  a  defined  file  of  the  external  tape. 
These  shock  data  can  be  read  back  with  User  Defined  Key  No.  2. 

A  numerical  output  of  all  analysis  data  can  be  retrieved  by 
pushing  User  Defined  Key  No.  3.  As  for  all  other  output  procedures 
the  primary  address  of  the  selected  output  device  (screen,  printer, 
or  plotter)  has  to  be  specified  here. 

User  Defined  Key  No.  4  makes  the  computer  load  part  II  of  the 
program  while  keeping  all  values  of  all  parameters  and  variables  In 
Its  memory. 

5.4.3.  Third  Itow  (No.  16  through  20) 

The  third  row  of  user  defined  keys  produces  a  linear  plot  of 
the  shock  sgectra.  The  amplitude  scale  is  automatically  adjusted  to 
the  output.  Any  range  of  analysis  frequencies  can  be  selected. 

User  Defined  Key  No.  16  displays  a  plot  of  pseudo-velocity 
vs.  analysis  frequency  (maximax) 

User  Defined  Key  No.  17  displays  a  plot  of  rel.  velocity 
vs.  analysis  frequency  (Residual) 

User  Defined  Key  No.  18  displays  a  plot  of  acceleration  vs. 
analysis  frequency. 

User  Defined  Key  No.  19  displays  a  plot  of  rel. 
displacements  vs.  analysis  frequency. 

User  Defined  Key  No.  20  displays  a  plot  of  phase  (0  to  180 
degree)  vs.  analysis  frequency. 

5.4.4.  Fourth  Row  (No.  6  through  10) 

A  logarithmic  plot  of  analysis  data  can  be  retrieved  by 
pushing  any  of  the  user  defined  keys  in  the  fourth  row.  User 
Defined  Key  No.  6  displays  the  Maximax  Shock  Spectrum,  No.  7  the 
Residual  Shock  Spectrum,  and  No.  8  both  of  them.  User  Defined  Key 
No.  9  allows  some  analysis  work  on  the  screen  in  determining  modal 
frequencies  and  appropriate  shock  spectrum  values  (the  Design  Shock 
Spectrum)  by  placing  the  two  crosshairs  at  points  of  interest.  The 
thumbwheels  on  the  keyboard  control  the  movement  of  the  crosshairs 
When  a  keyboard  key  is  pressed,  the  location  of  the  crosshair  is 
marked  and  the  appropriate  analysis  frequency  and  shock  spectrum 
value  is  printed  beside  the  graph-plot.  When  the  key  "0"  is 


* 


except  phase  is  always  displayed  between  0  and  180  degrees. 
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pressed,  the  crosshair  disappears  from  the  screen.  (Note:  this 
procedure  works  with  the  screen  only!) 


User  Defined  Key  No.  10  gives  a  plot  of  a  number  of  different 
shock  spectra  stored  on  the  external  tape.  Fach  spectrun  is  drawn 
with  a  different  dark  line  and  the  description  of  the  spectrum  is 
printed  wherever  the  crosshairs  have  been  placed. 

6.  Examples  of  Analysis  System  Output. 

6.1.  Figures  3  and  4  show  the  various  output  formats  which  can  be 
selected.  These  were  made  using  the  hard  copy  unit,  which  produces  a 
reduced-size  version  of  the  CRT  display.  Higher  quality  graphs  can  be 
obtained  from  the  plotter. 

6.2.  For  Fig.  3,  the  input  record  was  a  half-sine  acceleration  pulse. 

The  original,  integral  and  double  integral  are  shown  in  Fig.  3A  as  scaled 
acceleration,  velocity  and  displacement.  The  maximax  shock  spectrum  is 
shown  in  linear  plot,  without  grid-lines,  in  Fig.  3R,  and  in  log-log 
plot,  with  grid-lines,  in  Fig.  3C.  The  maximax  and  residual  are  shown  in 
log-log  plot,  with  factor-of-ten  grid-lines,  in  Fig.  3D.  Finally,  a 
tabular  listing  of  all  computed  data  is  shown  in  Fig.  3E. 

6.3.  The  input  record  for  Fig.  4  was  an  acceleration-time  history 
recorded  on  a  resilient  structure  aboard  the  Floating  Shock  Platform  with 
shot  standoff  of  20  ft.  The  acceleration,  velocity  and  displacement 
waveforms  are  shown  in  Fig.  4A,  utilizing  the  offset  correction  routine 
described  above:  i.e.,  requiring  the  average  velocity  near  the  400 
millisec.  mark  to  be  zero.  The  maximax  and  residual  spectra  are  given  in 
Fig.  4B,  and  again  in  Fig.  4C,  which  demonstrates  the  extraction  of 
design  shock  spectral  parameters.  Modal  frequencies  and  shock  spectrum 
values  are  lined  up  with  the  CRT  crosshairs,  and  their  values  tabulated. 
The  values  read  are  also  identified  on  the  graph.  The  tabular  listing  of 
computed  data  is  shown  in  Fig.  4D. 
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Fig.  3  —  Output  for  analysis  of  a  half-sine  acceleration  pulse 
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(d)  Maximax  and  reaidual  shock  spectra,  log-plot,  with  factor-of-ten  grid 
Fig.  3  —  Output  for  analysis  of  a  half-sine  acceleration  pulse 
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(e)  Tabular  listing  of  all  computed  data 
Fig.  3  —  Output  for  analysis  of  a  half-sine  acceleration  pulse 


.S3, 


(a)  Acceleration,  velocity,  and  diaplacement  time  histories 
Fia.  4  ~  Output  for  analysis  of  FSP  acceleration  record 


(2«ro  Ooaoirg)  Noturol  Fraquartcy  Hz 

(b)  Maximax  and  residual  shock  spectra 
Fig.  4  —  Output  for  analysis  of  FSP  acceleration  record 
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(d)  Tabular  Hating  of  all  computed  data 
Fig.  4  —  Output  fur  analysis  of  FSP  acceleration  record 
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Appendix  1 

List  of  Program  Variables 


1  A  Arroy  for  rim*  record  OIf1(N0J 
2:A2  Zero  offset  -  dotrjx  010(7.21 

3  B  Array  010(21 

4  B3  Vorioble  used  in  subroutine  TlOE 
5.B4  Vorioble  used  m  subroutine  TIOE 
6  C  Arroy  OIO(C01 

7 ; C0  Actual  number  of  analysis  frequencies  (defoult  901 
BiCI  Element  of  P'5  (first  analysis  frequency  per  frequency  rongel 
9 1 C2  Element  of  F5  (lost  onolysis  frequency  per  frequency  rongel 
10  C3  Oox  possible  analysis  frequencies  (depends  on  memory  ovailoble) 

110  I/O  device  oddness  (default  321 

1201  -ile  number  on  ext.  tope  drive  to  store  doto 

;3iD3  Pile  number  on  ext.  tope  drive  to  reod  doto  from 

14. E  Instantaneous  factor  for  sample  step 

15. EI  ."octor  for  somple  step  sitOin  frequency  ronge  lll-10Hzl 

16  E2  Poet  or  for  somple  step  within  frequency  ronge  [I(10-100hz> 

17. E3  Poctor  for  sample  step  within  frequency  ronge  1 1 ! ( 1 00- 1 000Hz l 

t8.F  Instantaneous  value  of  onolysis  frequency 

IQrFl  Lowest  onolysis  frequency  (default  II 
20:F2  Highest  onolysis  frequency  (default  10001 

21  ;F3  -owe*"  onolysis  frequency  within  one  of  the  3  frequency  ronges 

22;F4  Upper  onolysis  frequency  within  one  of  the  3  frequency  ronges 

23!F5  Frequency  arroy  O1M(C01 

24:F6  Lower  frequency  of  lin,  graph  (input  I 

25iF7  Upper  frequency  of  lin.  groph  (input! 

2SiF8  Cutoff  frequency  of  lowposs  filter 

27  G  Instontoneous  volue  of  occeloroiion  ( f2*pi»F»  f2ltXl 

28  G5  Acceleration  arroy  OIf1(C01 

29  H  Time  increment  of  digitized  time  record 

30  HI  Exponent 

3’  HZ  Time  increment  for  computation 

32  I  Counter 

33  1 1  Counter  to  find  Cl 

34- 12  Indicates  f requency  range  (I  to  II II 
35:13  Sample  number  when  V2  occurs 
36.14  Sample  number  when  VI  occurs 
37  K  Target  or^oy  for  groph  disp  loy 

38-M  Factor  of  F8  to  calculate  somple  frequency  withm  frequency  ronqe  I  (d 
e^oul >21 

39  A2  Poctor  of  F8  to  calculate  somple  frequency  withm  frequency  range  II  ! 
def oul t  21 

40  A3  Factor  of  F8  to  colculote  sample  frequency  wiih.n  frequency  range  11 ■ 
(def  oul 1  1 01 

4 1  L  Cali  brat  ion  level  in  i npu  t  un i t  s 
42 .  r*  Number  of  doto  module 

43  ni  tloximum  of  orroy  N  for  oufoscole 

44  HQ  Colibrotion  Hotrix  0Ih(7,21 

45  N0  Number  of  doto  points  on  record  (defoult  40961 
46:N  Number  of  doto  points  considered  during  colculot lon 

47  Nl  _ower  vert, col  level  of  log  groph  disploy 

48  N2  Upper  vert. col  level  of  log  groph  disploy 

49. N4  Number  of  doto  points  at  the  beginning  of  the  residuol  response 
50:N5  Number  of  doto  points  after  irunctot.on 
5’  N6  Lower  frequency  of  log  graph  display 
52  N7  Uoper  f requency  of  log  groph  display 

53. C  Indicotes  'he  port  of  the  program  currently  looded  (i  or  2i 
54  P  Selects  g- . d  on  log.  display  (P>!0  no  grid,  defoult  1 01 
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55. PI  Error  indicoror  if  >0  902 

5S  Q  Nu*!*ber  of  periods  considered  per  analysis  frequency 

57  S  Selecis  riAXIMAX-  and/or  RESIDUAL  SHOCK  SPECTRun ( Keys  6  to  9l 

58  R1  Selects  div.  linear  spectre  (user  defined  Keys  16  to  20) 

53  R2  Ind  cates  di splocenenT ( 1 ) , veloci ty 12)  or  accelerot ion(3)  record 

60  R3  Equcls  X3  before  N4  (initiol) 

61  .  R4  Equcls  X3  of  ter  NA  (residuall 

62  R5  Flag  for  of f set -procedure 

63  R6  Flog  for  deleting  calibration  routine 

64  R8  Flog  m  calbrotion  routine 

65  S  Calibration  factor 

66  SI  Calibrotion  factor  after  d i f f eren t i of i on  ond/or  integrotion 
67:S2  Calibration  foctor  of  ter  E-procedure 

63  T  I  ns t ont aneous  value  of  time  of  maximum  IT3  or  ^4) 

69 •T I  Beginning  of  time  record 

70  T3  Time  of  max  respo.nse  before  N4  'initioJ) 

7'  t4  Time  o‘  max.  response  after  N4  (residual) 

72  T5  Time  arroy  (contO'ns  max  volues  of  either  T3  or  T41  OlPltCB) 

73  U  Time  difference  m  oecimol  volue 

74  V  Instontoneous  volue  of  pseudvel oc i t y  (I2»pi*F>*X) 

75 : V '  Maximum  of  orroy  A 

76  V2  Mini  mum  of  array  A 

77  :V5  Veloci  *y  arroy  Oin(C0'i 
78 ; v  Foe  t  or 

79  W  Number  of  spectra  pro  plot  (Multi  Disploy  Mode)  and  flag  for  ws-aisplo 
y 

80  W2  Array.  DIM  W2(W1 ) 

81 iX  X-coordinate  of  pointer 

82  X  Instontoneous  value  of  displocement  tmoximum  of  R3  or  R4  i 

83:X‘  Deflection  ot  -0.54H2 

84:X2  Deflection  at  ♦0.54H2 

65; x3  Mox  response  of  analysis  frequency  F 

86;X5  Displacement  orroy 

87, Y  '' -coord  i  no  te  of  pointer 

88  Y  ■ ns f on roneous  volue  of  of  mox.  response  velocity  after  N4 
8S;Y5  Veloc.ty  orroy  DIM(C0) 

90. Z  Instontoneous  value  of  of  phase 
SliZS  Phase  orroy  D[M(C0) 

92: A*  String  for  groph  disploy 

93i3*  System  time  when  process-procedure  is  completed 
94 ! C*  Time  siring 

95.0*  Siring  to  address  do?a  module 

96  G*  String  for  gi-oph  disploy 

97: H*  String  for  graph  display 

98  L*  String  iq  checK  cond’ion  of  dato  module 

99 ■ ns  String  for  groph  disploy 

1 00  N*  Siring  for  groph  di^soloy 

101  P*  S*ring  for  groph  display 

102  S’ring  [  f  '0"  stops  anolysis  procedure) 

03 ;R*  "NRL  -  CODE  5842" 

104:T*  Contains  system  time 
05:V*  Contains  text  for  vert  level 
’  06  W*  Tiile  10  identify  record 
|07:X*  Segment  string 
'08  YS  Segment  sir  mg 
109  Z*  Segment  string 
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Appendix  2 
Program  Listing 

1  REM  ////////////////////////  eSPRGM/PART  I  \\\\\\\\\\\\\\\\\\\\\\\\\\\ 

2  INIT 

3  GO  TO  260 

4  IF  0=1  THEN  92 

5  GO  TO  7450 

8  IF  0=1  THEN  92 

9  GO  TO  7530 

12  IF  0=1  THEN  92 

13  GO  TO  8980 
16  GO  TO  190 
20  GO  TO  100 

24  IF  0=1  THEN  92 

25  R»1 

26  GO  TO  7690 

28  IF  0=1  THEN  92 

29  R»2 

30  GO  TO  7690 

32  IF  0=1  THEN  92 

33  R=3 

34  GO  TO  7690 

36  IF  0=1  THEN  92 

37  GO  TO  9510 

40  IF  0=1  THEN  92 

41  GO  TO  9780 

44  IF  0=2  THEN  90 

45  CHARSIZE  3 

46  IF  R6=0  THEN  94 

47  IF  R6*l  THEN  4200 

48  IF  0=2  THEN  90 

49  GO  TO  530 

50  DELETE  A 

52  IF  0=2  THEN  90 

53  GO  TO  2510 

54  DELETE  A 

56  IF  0=2  THEN  90 

57  GO  TO  5250 

60  IF  0=2  THEN  90 

61  R5=l 

62  CHARSIZE  3 

63  GO  TO  2410 

64  IF  0=1  THEN  92 

65  Rl=l 

66  GO  TO  6010 

68  IF  0=1  THEN  92 

69  Rl=4 

70  GO  TO  6010 

72  IF  0=1  THEN  92 

73  Rl=2 

74  GO  TO  6010 

76  IF  0=1  THEN  92 

77  Rl=3 

78  GO  TO  6010 

80  IF  0=1  THEN  92 

81  Rl=5 

82  GO  TO  6010 

90  PRINT  "You  got  to  PRESS  KEY  5  first !G" 

91  RETURN 

92  PRINT  "You  got  to  PRESS  KEY  4  firstlG" 

93  RETURN 

94  PRI  "LJJTHE  CALIBRATION  PROCEDURE  HAS  BEEN  DELETED  IN  ORDER  TO  SAUE"; 
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'?5  PRINT  '■  MEMORY.  IF  YOU  WANT  TO  REDEFINE  THE  CALIBRATION  MATRIX  M9 

96  PRINT  "RELOAD  THE  PROGRAM  (PRESS  USER  DEFINED  KEY  4  AND  5  IN  A  ROW" 

97  PRINT  "OR  PRESS  > AUTOLOADO I GG " 

99  return 

100  REM  »**»»*■»*•»»**»**  LOAD  ROUTINE 
110  IF  0=1  THEN  160 

120  DELETE  F5 , X5 , US , G5 , T5 , Y5 . Z5 , C , K , A 
130  R6=l 
140  0=1 

150  APPEND  "eSPRGM/PARTI"; 1 

160  PRINT  "LGSPRGM/PART  I  IS  LOADED!" 

170  DELETE  5380,20000 
130  RETURN 
190  IF  0=2  THEN  240 
200  0=2 

210  DELETE  270.20000 

220  APPEND  "eSPRGM/PARTII" j260 

230  RENUMBER  6000,10,260 

240  PRINT  "LGSPRGM/PART  II  IS  LOADED!" 

250  RETURN 

260  REM  x*»«»****»»****  DEFAULT  CONDITION  ****»»»»«■»»»*•»*»»»*»»«»»*»»»»* 
270  U$="TEST" 

280  0=1 
290  U1=0 
300  R5=0 
310  R6=l 
320  D=32 
330  Fl=l 
340  F2=1000 
350  F8=1000 
360  Nl=0.01 
370  N2=100 
380  N6=l 
390  N7=1000 
400  C0=90 
410  N0=4096 
420  Kl=2 
430  K2=2 
440  K3=10 
450  P=10 

460  DIM  M9(7,2) ,C(C0 ) ,A2<7,2) .R*( 15) ,T»( 18) ,B*( 18) ,C»( 18) 

470  i;IM  D*(8) ,H*(20) ,y$(22) ,G$<7) ,M*(7) ,N$(8) ,X$(22) ,P»(35) ,L*(4) ,Q$( 1 ) 
480  R4="NRL  -  CODE  5842" 

490  M9- 0 
500  A2=0 

510  PRINT  "  GSPRGM/PART  I  IS  LOADED  (DEFAULT  CONDITIONS)!" 

520  RETURN 

530  REM  »»**»**»*»■»*■»»•»»**  READ  RECORD 
540  DELETE  4200,5240 
550  R6=0 

560  PRINT  SD, 17: 1 .2, 1 .3 
570  CHARSIZE  3 

580  PRI  "LTHE  CURRENT  TITLE  IS  J  ENTER  A  NEW  TITLE  IF  NECESSA"; 

590  PRINT  "RY,  OTHERWISE  PRESS  'RETURN'. JC" 

600  INPUT  A4 

610  IF  A4=""  THEN  630 

620  W<V=AS 

630  DIM  A(N0) 

640  CHARSIZE  3 
650  IF  R5=0  THEN  660 
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660  U3«130 

670  IF  D<  >4  THEN  690 
680  U3«150 

690  PRINT  "JENTER  NUMBER  OF  DATA  MODULE  TO  READ  FROM:G"; 

700  INPUT  M 

710  IF  R5»l  THEN  740 

720  IF  M9<M,2)=0  THEN  2390 

730  S=M9(M,2) 

740  PRINT  86:"P" 

750  GO  TO  M  OF  760,780,800,820,840,860,880 
760  D*="I01B4096'' 

770  GO  TO  890 
780  D*=“I02B4096" 

790  GO  TO  890 
800  D«=''I03B4096" 

810  GO  TO  890 
820  D*="I04B4096" 

830  GO  TO  890 
840  D4="I05B4096" 

850  GO  TO  890 
860  D«="I06B4096" 

870  GO  TO  890 
880  D*="I07B4096” 

890  PRINT  e6;D« 

900  INPUT  e6!A 
910  P*="D" 

920  PRINT  e6!P4 
930  INPUT  86!?* 

940  A=»A-A2<M,2) 

950  H*-‘'SA" 

960  PRINT  BBiH* 

970  INPUT  86 iH* 

980  H1=VAL<H*) 

990  IF  HU12  THEN  1180' 

lOOO  GO  TO  HI  OF  1,1,1,1010,1030,1050,1070,1090,1110,1130,1150,1170 

1010  H»0.005 

1020  GO  TO  1200 

1030  H=0.ai 

1040  GO  TO  1200 

1050  H=0.025 

1060  GO  TO  1200 

1070  H=0.05 

1080  GO  TO  1200 

1090  H=0 . 1 

1100  GO  TO  1200 

1110  H=0.25 

1120  GO  TO  1200 

1130  H«0.5 

1140  GO  TO  1200 

1150  H»1 

1160  GO  TO  1200 

1170  H-2.5 

1180  PRINT  "ENTER  TIME  INCREMENT  IN  MS 

1190  INPUT  H 

1200  IF  R5=l  THEN  2430 

1210  CHARSIZE  2 

1220  CALL  ''TIME",T* 

1230  PRINT  8D: "LK";?* 

1240  PRINT  8D:T* 

1250  UINDOU  0,U3, 0,100 
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1260  yiEUPORT  0,U3, 0,100 
1270  IF  W1 >0  THEN  1350 
12B0  PRINT  eD,17!2,3 
1290  CHARSIZE  4 
1300  MOVE  8D:W3/2,98 
1310  FOR  1=1  TO  LEN<U*)/2 
1320  PRINT  @0: "H" ; 

1330  NEXT  I 
1340  PRINT  eD;W* 

1350  PRINT  eD,17!l .2,1 .8 

1360  CHARSIZE  2 

1370  Tl=H»AEiS(VAL<P*) ) 

1300  GO  TO  R2  OF  1390,1520,1660 

1390  GOSUEi  1840 

1400  S1=S 

1410  GOSUB  2020 

1420  PRINT  @D 

1430  CALL  "DIF3",A,A 

1440  S1=S1/H 

1450  GOSUB  1900 

1460  GOSUB  2020 

1470  CALL  '■DIF3",A,A 

1480  S1=S1  *l0''3/(H»9.a044) 

1490  GOSUB  i960 
1500  GOSUB  2020 
1510  GO  TO  1780 
1520  GOSUB  1900 
1530  S1=S 
1540  GOSUB  2020 
1550  PRINT 

1560  CALL  ••D1F3%A,A 
1570  S1=S1»10*3/<9,8044*H) 

1580  GOSUB  1960 

1590  GOSUB  2020 

1600  CALL  "INT  ”,A,A 

1610  CALL  "INT  ",A,A 

1620  S1=S1»1 . OE-3»H*2»9.8044 

1630  GOSUB  1840 

1640  GOSUB  2020 

1650  GO  TO  1780 

1660  GOSUB  1960 

1670  S1=S 

1680  GOSUB  2020 

1690  PRINT  8D: 

1700  CALL  "INT”, A, A 

1710  S1=S1»1 .0E-3»H»9.8044 

1720  GOSUB  1900 

1730  GOSUB  2020 

1740  CALL  "INT", A, A 

1750  S1=S1»H 

1760  GOSUB  1840 

1770  GOSUB  2020 

1780  VIEWPORT  0,W3, 0,100 

1790  WINDOW  0,W3, 0,100 

1800  MOVE  eD;0,0 

1810  PRINT  @D;  "*  oriq.  Tims  record 

1820  delete  a 
1830  RETURN 

1840  VIEWPORT  5,W3-30,4,32 
1850  V*="DISPLACEMENT" 


1860  G*="D  Max  :  " 

U370  M»="D  Min 
1880  N*="  [mm]" 

1390  RETURN 

1900  VIEWPORT  5,W3-30,36,64 
1910  V»="VELOCITY" 

1920  G*="V  Max  ; " 

1930  H$="V  Min 
1940  N«="  [M/sec]" 

1950  RETURN 

1960  VIEWPORT  5 , W3-30 , 68 , 96 
1970  V*="ACCELERATION" 

1980  Gi="A  Max  : " 

1990  M»="A  Min  i " 

2000  N$="  tgl" 

2010  RETURN 

2020  CALL  "MAX", A, VI ,14 
2030  CALL  "MIN", A, V2, 13 
2040  WINDOW  0,N0,V2,V1 
2050  IF  D=32  THEN  2110 
2060  MOVE  88:0,0 
2070  FOR  1=1  TO  NO 
2080  DRAW  eD:I,A<I) 

2090  NEXT  I 

2100  GO  TO  2120 

2110  CALL  "DISP",A 

2120  AXIS  8D:400,0 

2130  IF  R5=l  THEN  2340 

2140  FOR  1*0  TO  NO  STEP  2000 

2150  MOVE  8D: 1,0 

2160  PRINT  SD: "HJ",I»H; 

2170  NEXT  I 

2180  PRINT  88: "HHHHHHHHHHHtMS]" 

2190  MOVE  8D:4200,vT 

2200  PRINT  88:  USING  2240 : “J" ,G<,V1 »S1 ,N$i 
2210  MOVE  88:4200, VI 
2220  PRINT  88:  "J"; 

2230  PRINT  88:  USING  2240 : "J" ,M»,V2*S1 ,N$; 

2240  IMAGE  A . 7A , 68 , 1 8 . 1 1 A 
2250  MOVE  88: 0 ,0 .5*(V2+V1 ) 

2260  PRINT  88: "HH"  j 
2270  FOR  1*1  TO  LEN<V*)/2 
2280  PRINT  88:  "K"; 

2290  NEXT  I 

2300  FOR  1=1  TO  LEN<V$) 

2310  X«*SEG(V$, 1 , 1  ) 

2320  PRINT  eD:X*;"HJ"; 

2330  NEXT  I 
2340  RETURN 

2350  PRINT  8D:"JDATA  MODULE  ";M)"  IS  IN  GROUND  MODE  !!GG" 
2360  RETURN 

2370  PRINT  "JDATA  MODULE  ";Mj"  IS  IN  CALIBRATION  MODE  !!GG" 
2380  RETURN 

2390  PRINT  "JREAD  THE  CALIBRATION  STEP  FIRST !C" 

2400  RETURN 

2410  PRINT  "LOFFSET  -  CORRECTION:" 

2420  GO  TO  630 
2430  VIEWPORT  0,110,10,90 
2440  CALL  "INT",A,A 
2450  GOSUB  2020 
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2460  POINTER  X,Y,a* 

2470  A2(M,1)=M 

2480  A2<M,2)=Y/X+A2<M,2) 

2490  R5*0 
2S00  RETURN 

2510  REM  »**»*»*■»»»*»*»«*»•)*»*»  PROCESS  •»»*»»»«»•»»■*•»*»»»»*»»*»**«»«***»»»» 
2520  CHARSIZE  3 
2530  CALL  "TIME'',B« 

2540  PRINT  eOi  “L";B»}'‘J" 

2550  PRINT  @32,26:2 
2560  IF  F1<1  THEN  2620 
2570  IF  FI >1000  THEN  2620 
2580  IF  F2<1  THEN  2630 
2590  IF  F2>1000  THEN  2630 
2600  IF  F2<=F1  THEN  2630 
261 «  GO  TO  2660 
2620  IF  C0»1  THEN  2660 

2630  PRI  “ANALYSIS  FREQUENCY  IS  OUT  OF  RANGE  ENTER  NEW  VALUES  OF  F1,F2G" 
2640  INPUT  F1,F2 
2650  GO  TO  2530 

2660  DELETE  F5 , X5 , V5 ,GS , T5 , Y5 ,Z5 ,C , K , A 
2670  DIM  A(N0) 

2680  C3*INT( <MEMaRY/8-18)/8) 

2690  DIM  F5(CO),X5(CO>,V5<CO>,GS<CO),T5<CO>,Y5(CO>,Z5(CO),C(CO) 

2700  El=INT<1000/<Kl*Fa»H) ) 

2710  E2*INT<1000/<K2»F8*H) ) 

2720  E3=INT(1000/(K3»F8«H) ) 

2730  F5<1)=F1 

2740  IF  C0=1  THEN  2780 

2750  FOR  1=2  TO  CO 

2760  F5<I)=<F2/F1 )^<1/<C0-1 ) )«F5(I-1) 

2770  NEXT  I 
2780  X5=0.01 
2790  V5=0.0l 
2800  G5=0.01 
2810  T5=0,01 
2820  Y5=0.01 
2830  Z5=0.01 
2840  11=1 
2850  12=0 

2860  IF  FI <10  THEN  2900 
2870  12=1 

2880  IF  Fl<100  THEN  2900 

2890  12*2 

2900  GO  TO  3070 

2910  E=E1 

2920  GOSUB  3160 

2930  IF  II >C0  THEN  3010 

2940  GO  TO  3070 

2950  E=E2 

2960  GOSUB  3160 

2970  IF  I1>C0  THEN  3010 

2980  GO  TO  3070 

2990  E*E3 

3000  GOSUB  3160 

3010  DELETE  A 

3020  CALL  "TIME",T$ 

3030  COSUB  4030 
3040  PRINT  @32,26:0 

3050  PRINT  »0:"_MAX.  POSSIBLE  ANALYSIS  FREQUENCIES:  “jC3 
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3060  RETURN 
3070  F3=F5(11) 

3080  C1=I1 

3090  12=12+1 

3100  FOR  11=1  TO  CO 

3110  IF  F5<I1)>10^I2  THEN  3130 

3120  NEXT  II 

3130  F4=F5(I1-1) 

3140  GO  TO  12  OF  2910,2950,2990 

3150  RETURN 

3160  PRINT  86: "P";D$ 

3170  INPUT  86: A 
3180  A=A-A2(M,2) 

3190  IF  E<=1  THEN  3260 
3200  J=1 

3210  FOR  1=1  TO  NO  STEP  E 
3220  A<J)=A(I> 

3230  J=J+1 
3240  NEXT  I 
3250  GO  TO  3270 
3260  E=1 
3270  N5=N0/E 
3280  H2=H»E 

3290  GO  TO  R2  OF  3300,3340,3370 

3300  CALL  "DIF3",A,A 

3310  CALL  '•DIF3'',A,A 

3320  S2=S 

3330  CO  TO  3380 

3340  CALL  "DIF3",A,A 

3350  S2»S»H2 

3360  GO  TO  3380 

3370  S2=S*10*-3*H2*2*9.8044 

3380  C2=I1-1 

3390  P1»0.0062831853»F3*H2 
3400  FOR  I*C1  TO  C2 
3410  C<I)=2»C0S<P1 ) 

3420  IF  C1=C2  THEN  3440 
3430  P1=P1»<F4/F3)*(1/(C2-C1 ) ) 

3440  NEXT  I 

3450  IF  Pl<0.902  THEN  3500 
3460  IF  143.566/H2>F2  THEN  3500 

3470  PRINT  "LGWARNING:  FREQUENCIES  ARE  TOO  HIGH  FOR  YOUR  TIME  INCREMENT" 
3480  PRINT  "  ERRORS  ARE  OUTSIDE  THE  RANGE  -7  TO  *3  PERCENT  FOR" 

3490  PRINT  "  FREQUENCIES  ABOUE  " ; 1 43 . 566/H2 ; "  HERTZ" 

3500  FOR  I=C1  TO  C2 
3510  F=F5<I) 

3520  IF  F<=1  THEN  3550 
3530  N=1000/<H2»5»LGT(F) ) 

3540  IF  N<N5  THEN  3560 

3550  N=N5 

3560  X1=0 

3570  X2»0 

3580  R3=0 

3590  R4-0 

3600  IF  F*0  THEN  3630 
3610  N4=N-l-500/<F*H2> 

3620  IF  N4>3  THEN  3640 
3630  N4=3 

3640  FOR  J=1  TO  N-1 
3650  X3=C(I)»X2-X1+A( J) 
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3660  IF  J<N4  THEN  3710 

3670  IF  ABS(X3><ABS<R4)  THEN  3740 

3680  R4=X3 

3690  T4=J 

3700  GO  TO  3740 

3710  IF  ABS<X3)<ABS(R3)  THEN  3740 

3720  R3=X3 

3730  T3=J 

3740  X1=X2 

3750  X2=X3 

3760  NEXT  J 

3770  IF  ABS(R3) )ABS<R4)  THEN  3800 

3780  R3=R4 

3790  T3=T4 

3800  X=S2*R3 

3810  y=0 , 00628318»F»X 

3820  G=6. 28831 8»F»V/9. 8044 

3830  T=Tl+(T3-0 .5)*H2 

3840  Y=0  .  00628318*F)fS2*R4 

3850  Z=1 .OE-3»<T1+<T4-0.5)*H2)*F+0.25 

3860  Z=2-INT<Z) 

3870  IF  Y>0  THEN  3890 
3880  Z=2-0.5 

3890  IF  Z<=0.5  THEN  3910 

3900  Z=Z-1 

3910  Z=360*Z 

3920  F5<I)=F 

3930  X5<I)=X 

3940  y5(I)=y 

3950  G5<I)=G 

3960  T5(I)=T 

3970  Y5(I)=Y 

3980  Z5(I)=Z 

3990  PRINT  SDi  USING  4000 ; I ,F ,V,N; 

4000  IMAGE  20T,3D,4X,4D, 1D,4X,5D,2D,3X,6D 
4010  NEXT  I 
4020  RETURN 

4030  REM  »**»»»  SUBROUTINE  TIME  **»»»»»»'»»*x**»*»*»*»»**«»***innt*»»»»» 
4040  DIM  B(2) 

4050  C»=B* 

4060  FOR  1=1  TO  2 
4070  X*=SEG<C*, 1 1 ,8) 

4080  Y«=SEG<X*,4,2) 

4090  Z$=SEG(X«,7,2) 

410  0  B(  I  )=UAL(X*)-*-(VAL(  Y*)  »60-*-VAL<Z«>  )/360  0 
4110  C$=TS 
4120  NEXT  I 
4130  U*B(2)-B(1) 

4140  B3=<U-INT<U) )*60 
4150  B4=(B3-INT(B3))»60 

4160  PRINT  BD;  ■•__CALCULATION  TIME  (hr  ;  Min  ;  sec  )  ! 

4170  PRINT  80:  USING  4180  :  INT(U)  ,  “  : ’’.INTIBS)  ,  "  :  "  ,  INT(B4) 

4180  IMAGE  2D, A, 2D, A, 2D 
4190  RETURN 

4200  REM  CALIBRATION  xx-xxxxxxxxxxxxxxvxxxx'xxx-xwxxxx****# 

4210  DIM  A<N0) 

4220  CHARSIZE  3 

4230  PRINT  "LARE  YOU  DEALING  WITH  J  A  DISPLACEMENT  RECORD  (MM)"; 

4240  PRINT  "  (ENTER  1)  A  UELOCITY  RECORD  (M/SEC)  (ENTER  2)  OR  A"; 

4250  PRINT  "  ACCELERATION  RECORD  (G)  (ENTER  3)5 
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4260  INPUT  R2 

4270  PRI  "  DO  YOU  WANT  TO  READ  A  CALIBRATION  STEP  (ENTER  1)  OR"; 
4280  PRINT  "  TO  ENTER  CALIBRATION  FACTOR  (ENTER  2)G 
4290  INPUT  R8 
4300  M9=0 

4310  IF  R8=2  THEN  4960 
4320  FOR  M=1  TO  7 

4330  PRI  "LENTER  CALIBRATION  LEUEL  FOR  CHANNEL  "  ; M ;  "  IN  INPUT  UNIT5G;"; 

4340  INPUT  L 

4350  GOSUB  4750 

4360  IF  S4=0  THEN  4390 

4370  IF  S4=7  THEN  4420 

4380  GO  TO  4450 

4390  GOSUB  2310 

4400  CALL  "WAIT", 2 

4410  GO  TO  4730 

4420  GOSUB  2330 

4430  CALL  "WAIT", 2 

4440  GO  TO  4730 

4450  PRINT  @6; "P" 

4460  GO  TO  M  OF  4470,4490,4510,4530,4550,4570,4590 
4470  D4='*I01B4096'‘ 

4480  GO  TO  4600 
4490  D*="I02B4096" 

4500  GO  TO  4600 
4510  D*="I03B4096" 

4520  GO  TO  4600 
4530  D«=''I04B4096“ 

4540  GO  TO  4600 
4550  0*="I05B4096" 

4560  GO  TO  4600 
4570  DS="I06B4096" 

4580  GO  TO  4600 
4590  D«="I07B4096" 

4600  PRINT  @6:0* 

4610  INPUT  e6:A 

4620  WINDOW  0,5000,-512,512 

4630  VIEWPORT  0,130,0,100 

4640  PRINT  "LCHANNEL  ";M;"  JJSET  CROSSHAIR  TO  UPPER  CALIBRATION  LEVELG" 
4650  CALL  "DiSP“,A 
4660  AXIS 

4670  POINTER  XI , Y1 ,0* 

4680  PRINT  ""J’JJ  AND  LOWER  CALIBRATION  LEVELG!" 

4690  POINTER  X2,Y2,Q* 

4700  IF  Y2>Y1  THEN  5230 

4710  M9<M,1)=M 

4720  M9(M,2)*L/(Y1-Y2) 

4730  NEXT  M 
4740  RETURN 

4750  GO  TO  M  OF  4760,4780,4800,4820,4840,4860,4880 
4760  L*="A01" 

4770  GO  TO  4890 
4780  L»="A02" 

4790  GO  TO  4890 
4800  L*="A03" 

4810  GO  TO  4890 
4820  L*="A04" 

4830  GO  TO  4890 
4840  L*="A05'’ 

4850  GO  TO  4890 
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4860  L$="A06" 

4870  GO  TO  4890 
4880  Li="A07" 

4890  PRINT  86; Li 
4900  INPUT  86;Li 
4910  Sl=INT(VAL<Li)/l 000) 

4920  S2=INT( (UAL<Li)-Sl *1 000 )/l 00 ) 

4930  S3=INT( (UAL (Li )-Sl *1  000-32*1 00 '/I  0  ) 

4940  S4=UAL(Li)-Sl *1 000-32*1 00-33*10 

4950  RETURN 

4960  PRINT  "LENTER." 

4970  FOR  M=1  TO  7 

4980  GOSUB  4750 

4990  IF  34=0  THEN  5020 

5000  IF  34=7  THEN  5040 

5010  GO  TO  5060 

5020  GOSUB  2310 

5030  GO  TO  5200 

5040  GOSUB  2330 

5050  GO  TO  5200 

5060  PRl  "JCALIBRATION  FACTOR  ( INPUT-UNITS/ 1 UOLT >  FOR  CHANNEL  "jMj"; 

5070  INPUT  L 

5080  GO  TO  34  OF  5090,5110,5130,5150,5170,5190 

5090  M9(M,2)=L/512 

5100  GO  TO  5200 

5110  M9(M,2)=L*0 .2/512 

5120  GO  TO  5200 

5130  H9(M,2)=L*5/512 

5140  GO  TO  5200 

5150  M9(M,2)=L*0 . 1/512 

5160  GO  TO  5200 

5170  M9(M,2)=L*2/512 

5180  GO  TO  5200 

5190  M9(M,2)=L*0 .5/512 

5200  M9(M,1)=M 

5210  NEXT  M 

5220  RETURN 

5230  PRINT  "'JJJJJ  UPPER  CALIBRATION  LEVEL  FIRST!'  TRY  IT  AGAINIIG" 

5240  GO  TO  4670 

5250  REM  »**»■»■**»*»***•»■*•»*■»*■**•*■»  WRITE  RECORD  »»*•»•»*•**•»**»»**•*«»»*****»* 
5260  DIM  A(N0) 

5270  CHARSIZE  3 

5280  PRINT  "LENTER  FILE  NUMBER  (DEVICE  NR. 3)  TO  WRITE  RECORD  ONG : 

5290  INPUT  D3 
5300  FIND  e3:D3 
5310  PRINT  86; "P";Di 
5320  INPUT  86;A 
5330  A=A-A2(M,2) 

5340  WRITE  83 ; Wi,R2,H,Tl ,M,M9(M,2) ,A2(M.2> ,A 
5350  PRINT  "JRECORD  IS  STORED  IN  FILE  ♦C'jDS 
5360  DELETE  A 
5370  RETURN 


6000  REM  SPRGM/PART  II  FILE  3 

6010  REM  ****■»«»»*»***■»*»»*»****  LIN.  SPECTRA  *«***»*»»■»**•»»»»*■»«******* 
6020  DELETE  K 
6030  DIM  K(C0) 

6040  Ai="MAXIMAX  SHOCK  SPECTRUM" 

6050  GO  TO  R1  OF  6060,6080,6100,6120,6150 
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6060  K=A6S(V5> 

6070  GO  TO  6170 
6080  K»ABS(GS) 

6090  GO  TO  6170 
6100  K»ABS(XS> 

6110  GO  TO  6170 

6120  A*=‘’RESIDUAL  SHOCK  SPECTRUM” 

6130  K^ABSCYS) 

6140  GO  TO  6170 
6150  A4= "PHASE" 

6160  K-Z5 

6170  CALL  "MAX”, K, Ml ,14 
6180  CHARSIZE  3 

6190  PRINT  "LJFREQUENCY  RANGEG:  ”; 

6200  INPUT  F6,F7 
6210  IF  D=4  THEN  6240 

6220  VIEWPORT  8»1 . 79 , 1 30-4*1 , 79 ,3»2 . 82 , 1 0 0-3*2 . 82 
6230  GO  TO  6260 

6240  VIEWPORT  8*1.79,150-4*1.79,3*2.82,100-3*2.82 

6250  PRINT  @0,17:1.2,1.8 

6260  CALL  "TIME”,T4 

6270  CHARSIZE  1 

6280  PRINT  @0: "L"}R$ 

6290  PRINT  eO'.T* 

6300  IF  Rl=5  THEN  6510 

6310  IF  ABS<M1)<10  THEN  6350 

6320  IF  ABS(M1)<100  THEN  6390 

6330  IF  ABS(M1)<1000  THEN  6430 

6340  IF  ABS(M1)>1000  THEN  6470 

6350  M1=«INT(M1  )  +  l 

6360  WINDOW  F6,F7,0,M1 

6370  AXIS  eD:<F7-F6)/10,0.5 

6380  GO  TO  6580 

6390  M1=INT(M1/10)*10+10 

6400  WINDOW  F6,F7,0,M1 

6410  AXIS  BD: (F7-F6)/10 ,5 

6420  GO  TO  6580 

6430  M1=INT(M1/100)*100+100 

6440  WINDOW  F6,F7,0,M1 

6450  AXIS  8D: (F7-F6)/10,50 

6460  GO  TO  6580 

6470  M1*INT<M1/1000)«1000+1000 

6480  WINDOW  F6,F7,0,M1 

6490  AXIS  SD: (F7-F6)/10,500 

6500  GO  TO  6580 

6510  WINDOW  F6,F7, 0,180 

6520  AXIS  SD: (F7-F6)/10,45 

6530  MOVE  eD;F6,180 

6540  DRAW  eD;F7,180 

6550  DRAW  eD;F7,0 

6560  MOVE  BD: (F7+F6)/2,ia0 

6570  CO  TO  6610 

6580  MOVE  SD:F6,M1 

6590  DRAW  ffD;F7,Ml 

6600  DRAW  SD:F7,0 

6610  PRINT  BD, 17:2,3 

6620  CHARSIZE  4 

6630  COSUB  6700 

6640  IF  WDO  THEN  6800 

6650  A«>W« 
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6660  PRINT  8D,17! 1 .2,1 .8 

6670  CHARSIZE  2 

6680  GOSUB  6700 

6690  GO  TO  6800 

6700  IF  Rl=5  THEN  6730 

6710  MOVE  90: <F7+F6)/2,M1 

6720  GO  rO  6740 

6730  MOVE  90: (F7+F6)/2, 180 

6740  PRINT  90: “KK " > 

6750  FOR  1=1  TO  LEN<A*)/2 
6760  PRINT  90: "H" ; 

6770  NEXT  I 
6780  PRINT  90:A$i 
6790  RETURN 

6800  PRINT  90,17:1 .2,2 
6810  CHARSIZE  2 

6820  H$="NaTural  Frequency  Hz" 

6830  MOVE  90: (F7+F6)/2,0 
6840  PRINT  90:"JJJ"j 
6850  FOR  1=1  TO  LEN(H*)/2 
6860  PRINT  90: "H"  j 
6870  NEXT  I 
6880  PRINT  90;H«; 

6890  GO  TO  R1  OF  6900,6940,6960,6920,6980 
6900  V$="Pseudowelocity  M/sec“ 

6910  GO  TO  7010 

6920  V$="Rel.  Velocity  M/sec" 

6930  GO  TO  7010 

6940  V$="Abs.  Acceleration  q" 

6950  GO  TO  7010 

6960  V«="Rel.  Displacewent  mm" 

6970  GO  TO  7010 
6980  V1i="  Phase  deg" 

6990  MOVE  9DiF6,90 

7000  GO  TO  7020 

7010  MOVE  9D:F6,0 .5*INT<M1  ) 

7020  PRINT  9D:"HHHHHH"; 

7030  FOR  1=1  TO  LEN<V*)/2 
7040  PRINT  9D:  "K  "  > 

7050  NEXT  I 

7060  FOR  1=1  TO  LEN<V*) 

7070  X$=SEG(V$,I,1  ) 

7080  PRINT  8D:X«; ”HJ"  j 
7090  NEXT  I 

7100  FOR  I=F6  TO  F7  STEP  <F7-F6)/10 
7110  MOVE  9D:I,0 
7120  PRINT  9D: "HJ" ; INT( I) j 
7130  NEXT  I 

7140  IF  Rl=5  THEN  7270 
7150  IF  ABS(M1)<=10  THEN  7190 
7160  IF  ABS(M1><=100  THEN  7210 
7170  IF  ABS(M1)<=1000  THEN  7230 
7180  IF  ABS<M1)>1000  THEN  7250 
7190  FOR  1=0  TO  ABS(Ml)  STEP  1 
7200  GO  TO  7280 

7210  FOR  1=0  TO  ABS(Ml)  STEP  10 
7220  GO  TO  7280 

7230  FOR  1=0  TO  ABS(Ml)  STEP  100 
7240  GO  TO  7280 

7250  FOR  1=0  TO  ABS<M1)  STEP  1000 


7260  GO  TO  7280 

7270  FOR  1=0  TO  180  STEP  45 

7280  MOVE  @D:F6,I 

7290  PRINT  80 ! "HHH" ;I ; 

7300  NEXT  1 
7310  MOVE  8D:F6,0 
7320  PRINT  80,17:1  ,2 
7330  A*=''K(Zero  OaMping)" 

7340  PRINT  80: 

7350  1=1 

7360  FOR  I=I  TO  CO 

7370  IF  F5(I)>0.01  THEN  7390 

7380  NEXT  I 

7390  MOVE  eD:F5<I) ,ABS<K<I) ) 

7400  FOR  I=I  TO  CO 

7410  IF  F5(I>=0.01  THEN  7360 

7420  DRAU  80 : F5< I ) , ABS<K < I)  ) 

7430  NEXT  I 
7440  RETURN 

7450  REM  URITE  OATA  tt******************************** > 

7460  CHARSIZE  3 

7470  PRINT  "JLGENTER  FILE  NUMBER  <OEVICE  NR. 3)  TO  WRITE  OATA  ON: 

7480  INPUT  01 
7490  FINO  83:01 

7500  URITE  83:C0,U*,F5,X5,V5,G5,T5,Y5,Z5 
7510  PRINT  "RESULTS  STOREO  IN  FILEG  ♦";01 
7520  RETURN 

7530  REM  xxxxxxxxxKxxxxxx**  REAO  OATA 
7540  CHARSIZE  3 

7550  PRINT  "JLGENTER  FILE  NUMBER  (OEVICE  NR. 3)  TO  REAO  BATA  FROM:  "j 
7560  INPUT  03 

7570  BELETE  F5 , X5 , V5 ,G5 , T5 , Y5 ,25 ,K  ,  A 
7580  FINO  83:03 
7590  READ  83:00 
7600  FINO  83:03 

7610  DIM  F5<C0) ,X5(C0) ,V5<C0) ,G5<C0) ,T5(C0) ,Y5(C0) ,Z5(C0) ,K(C0) 

7620  READ  83 : CO , WS , F5 , X5 , V5 , G5, T5 , Y5, Z5 

7630  IF  W1 >0  THEN  7650 

7640  PRINT  "DATA  ARE  IN  MEMORY'G* 

7650  RETURN 

7660  REM  LOG.  SPECTRA 

7670  VIEWPORT  5»1 . 79 , 1 1 0-4«l . 79 , 1 0+3*2 . 82 , 1 00-3*2 . 82 
7680  GO  TO  7740 

7690  VIEWPORT  8*1 . 79 , 1 30-4*1 . 79 , 3*2 . 82 , 1 00-3*2 . 82 
7700  IF  0=4  THEN  7720 
7710  GO  TO  7740 

7720  VIEWPORT  8*1.79,150-4*1.79,3*2.82,100-3*2.82 

7730  PRINT  80,17:1.2,1.8 

7740  CHARSIZE  1 

7750  DELETE  K 

7760  DIM  K<C0) 

7770  CALL  "TIME",T$ 

7780  PRINT  SO: "L"jR* 

7790  PRINT  BD:Ti> 

7800  REM 

7810  WINDOW  LGT(N6) ,LGT(N7) ,LGT(N1 ) ,LGT(N2> 

7820  IF  P)10  THEN  7980 
7830  H1=INT(LGT(N1 ) ) 

7840  FOR  I=10*H1  TO  10*(H1+1)  STEP  lO^HlxP 
7850  MOVE  8D:LCT<N6) ,LCT(1) 


7860  DRAW  8D : LGT ( N7 > ,LGT < I ) 

7870  NEXT  I 
7880 

7890  IF  H1<=LGT(N2)  THEN  7840 
7900  H1=LGT<N6) 

7910  FOR  I  =  10''H1  TO  10*<H1+1)  STEP  t0*Hl»P 
7920  MOVE  8D:LGT(I) ,LGT(N1 ) 

7930  DRAW  8D : LGT < I ) , LGT ( N2 ) 

7940  NEXT  I 
7950  H1=H1+1 

7960  IF  H1<=LGT<N7)  THEN  7910 

7970  GO  TO  8020 

7980  AXIS  8D; 1 ,1 ,0,-2 

7990  MOVE  8D:LGT<N6) ,LGT(N2) 

8000  DRAW  8D:LGT<N7) ,LGT(N2) 

8010  DRAW  8D:LGT<N7) ,LGT<N1 ) 

8020  PRINT  80,17:2, 3 
8030  CHARSIZE  4 

8040  GO  TO  R  OF  8050,8070,8090 
8050  A*="MAXIMAX  SHOCK  SPECTRUM" 

8060  GO  TO  8100 

8070  Ai="RESIDUAL  SHOCK  SPECTRUM- 
8080  GO  TO  8100 

8090  A$="MAXIMAX  and  RESIDUAL  SHOCK  SPECTRUM" 

8100  GOSUB  8170 

8110  IF  U1 >0  THEN  8240 

8120  A$=U« 

8130  PRINT  8D,17!l .2,1 .8 
8140  CHARSIZE  2 
8150  GOSUB  8170 
8160  GO  TO  8240 

8170  MOVE  @D:LGT<10"<0.5*LGT<N6»N7))),LGT(N2> 
8180  PRINT  8Di "KK") 

8190  FOR  1=1  TO  LEN<A«)/2 
8200  PRINT  8D: "H"; 

8210  NEXT  I 
8220  PRINT  SDjA*; 

8230  RETURN 

8240  GO  TO  R  OF  8250,8270,8290 
8250  V*="Ps*udo velocity  n/sec" 

8260  GO  TO  8300 

8270  V4="Rol.  Velocity  «/sec" 

8280  GO  TO  8300 

8290  V«=”Velocity  n/sec  " 

8300  MOVE  8D!LGT<N6),LGT<10"<0.5»LGT(N2*N1))) 

8310  PRINT  8D,17! 1 .5,2.5 

8320  CHARSIZE  2 

8330  PRINT  8D: "HHHHHHH" j 

8340  FOR  1=1  TO  LEN(V«)/2 

8350  PRINT  8D!  "K"; 

8360  NEXT  I 

8370  FOR  1=1  TO  LENCV*) 

8380  X«=SEG(V«,I,1 ) 

8390  PRINT  8D:Xi;"HJ"; 

8400  NEXT  I 

8410  MOVE  BD!LGT( 10*(0 .5»LCT<N6»N7) ) ) ,LGT(N1 ) 
8420  HO»"Natural  Frequency  Hz" 

8430  FOR  1=1  TO  LEN<H$)/2 
8440  PRINT  BD: "H"  ; 

8450  NEXT  I 
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3460  PRINT  SD: 

8470  H1=0 

8480  FOR  I=LGT(N1)  TO  LGT(N2)  STEP  lO^Hl 
8490  MOVE  eD!LGT(N6),I 

8500  PRINT  80:  "HHHHH"j"10*'';INT<Ht+LGT<Nl>) 
8510 

8520  NEXT  I 
8530  H1=:0 

8540  FOR  I=LGT(N6)  TO  LGT<N7)  STEP  10*H1 
8550  MOVE  80: I,LGT(N1 ) 

8560  PRINT  80: "HHJ";"10*“;INT(H1+LGT(N6>) 
8570  H1=H1+1 
8580  NEXT  I 

8590  MOVE  eO:LGT<N6) ,LGT<N1  ) 

8600  IF  Dc32  THEN  8620 
8610  PRINT  80,17:1,2 
8620  CHARSIZE  2 
8630  A$s"<Zero  Oanping)” 

8640  PRINT  80:“JJJ"jA* 

8650  GOSUB  8770 
8660  GOSUB  8820 
8670  IF  R=3  THEN  8690 
8680  RETURN 
8690  R=R+1 

8700  IF  004  THEN  8740 

8710  PRINT  "LYOU  MAY  CHANGE  THE  PEN  COLOR  G" 
8720  PRINT  •  Press  RETURN  to  continue" 

8730  INPUT  Ai 
8740  GOSUB  8770 
8750  GOSUB  8820 
8760  RETURN 

8770  GO  TO  R  OF  8780,8800,8780,8800 
8780  K=V5 
8790  GO  TO  8810 
8800  K=Y5 
8810  RETURN 
8820  1=1 

8830  FOR  I=I  TO  CO 
8840  IF  F5(I)>0.01  THEN  8870 
8850  NEXT  I 
8860  RETURN 

8870  MOVE  eD:LGT<F5(I) ) ,LCT<ABS<K<I)  )  ) 

8880  FOR  I=I  TO  CO 
8890  IF  F5<I)=0.01  THEN  8830 
8900  IF  R<*3  THEN  8920 
8910  DASH  170 

8920  DRAW  8D : LGT < F5< I ) ) , LGT( ABS<K < I ) )  ) 

8930  NEXT  I 

8940  IF  U1 >0  THEN  8960 
8950  DASH  0 
8960  HOME 
8970  RETURN 

8980  REM  »»»»»«»«»»*»»»  NUMERICAL  OUTPUT 
8990  WINDOW  0,130,0,100 
9000  VIEWPORT  0,130,0,100 
9010  11=0 
9020  12=1 
9030  N«1 

9040  IF  D=32  THEN  9070 
9050  IF  D=40  THEN  9130 


9060  PRINT  9D,17i 1 .5,2.5 
9070  CHARSIZE  3 
9080  PRINT  9D:“U" 

9090  HOVE  90:65,99 
9100  FOR  1=1  TO  LEN<U*)/a 
9110  PRINT  90: "H"; 

9120  NEXT  I 
9130  PRINT  80:U« 

9140  CHARSIZE  1 
9150  CALL  ”TIHE",T$ 

9160  PRINT  90:"*J“jR*,“  “jT*,“JI" 

9170  IF  0-40  THEN  9190 

9180  PRINT  90:  USING  9190 : "KKK “PAGE* ; 12; "JJ " 

9190  IHAGE  3A,12aX,4A,20,2A 
9200  CHARSIZE  2 

9210  PRINT  90:"  FREQ  SHOCK  SPECTRUM  FOURIER  TRANS* 

9220  PRINT  90:"  (HZ)  <MH)  (H/S>  (G)  (MS)  (M/S)  (DEC)_* 

9230  IMAGE  40 . 1 0, X , 50 . 20 , X , 40 . 20,X,5D . 10,X, 40 . 2D,X, 40 . 2D,2X, 40 
9240  IMAGE  65T , 40 . 1 D , X ,50 . 20 ,X , 40 . 20 ,X ,50 . 10 , X , 40 . 20 ,X , 40 . 20 ,2X , 40 
9250  FOR  I-l+Il  TO  CO 
9260  IF  I>50«N  THEN  9300 

9270  PRINT  90:  USING  9230 : F5( I ) ,X5( I ) , V5( I ) ,G5 ( I ) , T5 ( I ) , Y5( I ) , Z5( I ) 

9280  NEXT  I 
9290  RETURN 

9300  IF  0=40  THEN  9350 
9310  CHARSIZE  1 
9320  PRINT  90: “"JJJ21" 

9330  CHARSIZE  2 
9340  PRINT  90: *J" 

9350  FOR  I-1+50»*N  TO  CO 
9360  IF  I>(N-:-l)«50  THEN  9430 
9370  IF  0=40  THEN  9400 

9380  PRINT  90:  USING  9240 : F5( I ) , X5( I ) , V5( I ) ,G5( I ) , T5( I ) , Y5( I ) , Z5( I ) 

9390  GO  TO  9410 

9400  PRINT  90:  USING  9230 : F5 ( I ) , X5( I ) , V5( I ) ,G5( I ) ,T5( I ) , Y5( I ) , Z5 ( I ) 

9410  NEXT  I 

9420  RETURN 

9430  IF  0=40  THEN  9450 

9440  CALL  "WAIT", 4 

9450  11=100 

9460  I2-I2>1 

9470  N-N+2 

9480  IF  0=40  THEN  9500 
9490  CO  TO  9040 
9500  GO  TO  9250 

9510  REM  ANALYSIS  »»«*««*««««»«»•«»««««*«*««*«*««« 

9520  R=3 

9530  COSUB  7670 

9540  R-1 

9550  MOVE  90:3.11,2 

9560  PRINT  90;  USING  *2A, 5X ,4A,3X,7A* : "NR * , " t HZ  1  * , " CM/SEC) " 

9570  POINTER  X,Y,Q» 

9580  IF  a»-"0"  THEN  9750 

9590  MOVE  90:X-LGT(1 . 022) ,Y+LGT(1 . 022) 

9600  CHARSIZE  2 

9610  PRINT  90; "J*" jR 

9620  COSUB  9650 

9630  R=R4-1 

9640  CO  TO  9570 

9650  VIEWPORT  106,130,0,100 
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9660  UINDOU  0,32,0,100 
9670  HOVE  80:0,94 
9680  FOR  I»1  TO  R 
9690  PRINT  8D:"JJ''i 
9700  NEXT  I 

9710  PRINT  80:  USING  9720  :  •* J*  "  ,R ,  “  :  " ,  10*X ,  1  0*Y 
9720  IMAGE  2A , 20 , A , 2X , 40 . 1 0 , 3X ,30 . 0 

9730  VIEWPORT  5»1 . 79 , 11 0-4»l . 79 , 1 0+3»2 . 82 , 1 00-3»2 . 82 
9740  WINDOW  LGT(N6) ,LGT<N7) ,LGT<N1 ) ,LGT<N2) 

9750  DASH  0 
9760  HOME 
9770  RETURN 

9780  REM  »»»*»»«»»»«**»»»»»*»»  MULTI  DISPLAY  »* »»*«»« »*»»»*»»»»»»*»*»*»* 
9790  CHARSIZE  3 
9800  0-32 

9810  PRINT  "LJENTER  NUMBER  OF  DESIRED  SPECTRAG;  " ; 

9820  INPUT  wi 
9830  DIM  W2(W1) 

9840  PRINT  "  ENTER  FILE  NUMBERS  ON  DEVICE  NR . 3G : 

9850  INPUT  U2 

9860  PRINT  “_ENTER:  1  FOR  MAXIMAX-SPECTRA_2  FOR  RESIDUAL-SPECTRAg :  "> 

9870  INPUT  R 

9880  D3-W2<n 

9890  COSUB  7570 

9900  GQSUB  7690 

9910  COSUB  10060 

9920  IF  Wl-1  THEN  10050 

9930  FOR  10-2  TO  Wl 

9940  D3=W2<I0) 

9950  COSUB  7570 
9960  DASH  11«I0 
9970  GOSUB  8770 
9980  GOSUB  8820 
9990  GOSUB  10060 
10000  NEXT  10 
10010  DASH  0 
10020  Wl-0 

10030  DELETE  W2 , X , Y , 1 0 , D3 , Q» 

1  0040  W*-"TEST'' 

10050  RETURN 

10060  POINTER  X,Y,Q* 

10070  MOVE  8D;X,Y 
10080  DRAW  8D:X+0.25,Y 
10090  PRINT  BD: "  " jW9 
10100  RETURN 
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1  REM  /////////////////  eSPRGM/PART  III  \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\ 

2  INIT 

3  GO  TO  260 

4  IF  0»3  THEN  94 

5  GO  TO  7450 

8  IF  0=3  THEN  94 

9  GO  TO  7530 

12  IF  0=3  THEN  94 

13  GO  TO  8980 
16  GO  TO  190 
20  GO  TO  100 

24  IF  0=3  THEN  94 

25  R  =  1 

26  GO  TO  7690 

28  IF  0=3  THEN  94 

29  R-2 

30  GO  TO  7690 

32  IF  0=3  THEN  94 

33  R=3 

34  GO  TO  7690 

36  IF  0=3  THEN  94 

37  GO  TO  9510 

40  IF  0=3  THEN  94 

41  CO  TO  9780 
45  RETURN 

48  IF  0=2  THEN  90 

49  GO  TO  510 

50  DELETE  A 

52  IF  0=2  THEN  90 

53  GO  TO  1780 

54  DELETE  A 

64  IF  0=3  THEN  94 

65  Rl=l 

66  GO  TO  6010 

68  IF  0=3  THEN  94 

69  Rl=4 

70  GO  TO  6010 

72  IF  0=3  THEN  94 

73  Rl*2 

74  GO  TO  6010 

76  IF  0=3  THEN  94 

77  Rl=3 

78  GO  TO  6010 

80  IF  0=3  THEN  94 

81  Rl=5 

82  GO  TO  6010 

90  PRINT  "Too  got  to  PRESS  KEY  5  firstlC" 

92  RETURN 

94  PRINT  "You  got  to  PRESS  KEY  4  first !G* 

96  RETURN 

100  REMs»«*s*s***ssss*s  LOAD  ROUTINE 
no  IF  0=3  THEN  170 

120  DELETE  F5 , X5 , V5 , G5 , T5 , Y5 ,Z5,C,C0 ,K , A 
130  C0=90 
140  0=3 

150  APPEND  "9SPRGM/PARTIII“;1 

160  PRINT  "JCSPRGM/PART  III  IS  LOADED!" 

170  DELETE  5140,20000 
180  RETURN 
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190  IF  0-2  THEN  240 
200  0-2 

210  DELETE  270 >20000 

220  APPEND  "9SPRGM/PARTI1" j260 

230  RENUMBER  6000 >10 >260 

240  PRINT  "JGSPRGM/PART  II  IS  LOADED!" 

250  RETURN 

260  REM*»*«»«**««»**«  DEFAULT  CONDITION 
270  M*="TEST" 

280  0-3 
290  Ul-0 
300  D-32 
310  Fl-1 
320  F2-1000 
330  F8-1000 
340  Nl-0.01 
350  N2-100 
360  N6-1 
370  N7-1000 
380  CO-90 
390  NO-4096 
400  Kl-2 
410  K2-2 
420  K3-10 
430  P-10 

440  R«-"NRL  -  CODE  5842" 

450  DIM  M9<7>2) >A2<7>2) >C(C0) 

460  M9-0 
470  A2-0 

480  DIM  D0(8) >H$<20) >V*<22) >G$<7>  >M»<7>  >N«(8>  >X*(22) ,P*(35) >L*(4),Q*<1 ) 
490  PRINT  “  GSPRGM/PART  III  IS  LOADED  (DEFAULT  CONDITIONS)!" 

500  RETURN 

510  REM  READ  RECORD»»«***»4»***»«»*»««»**»'*'>(*’»«*(****4(*« 

520  CHARSIZE  3 

530  PRINT  "LJENTER  FILE  NUMBER  (ON  DEVICE  3)  TO  READ  RECORD  FROMGs  "} 

540  INPUT  dT 
550  DIM  A(N0) 

560  U3-130 

570  IF  D04  THEN  600 
580  U3-150 

590  PRINT  8D>17! 1 .2>1 .8 
600  CHARSIZE  2 
610  CALL  "TIME">T$ 

620  PRINT  6D: "UK">R» 

630  PRINT  eD:T« 

640  FIND  83;D1 

650  READ  83;U«>R2>H>T1 >M 

660  FIND  «3!D1 

670  READ  93:U*>R2>H>T1 >M>M9(M>2) >A2(M>2) >A 
680  M9(M>1)-M 
690  A2(M>1)-M 
700  S-M9(M>2> 

710  WINDOW  0>W3>0>100 
720  VIEWPORT  0>W3>0>100 
730  IF  W1>0  THEN  810 
740  PRINT  9D>17:2,3 
750  CHARSIZE  4 
760  MOVE  8D;W3/2>98 
770  FOR  I-l  TO  LEN(W0)/2 
780  PRINT  9D:"H"} 
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790  NEXT  I 
800  PRINT  eOiU* 

810  PRINT  @0,17:1.2,1.8 
820  CHARSIZE  2 

830  GO  TO  R2  OF  840,970,1110 

840  GOSUB  1280 

850  S1»S 

860  GOSUB  1460 

870  PRINT  eD:“»" 

880  CALL  "DIF3",A,A 

890  S1=S1/H 

900  GOSUB  1340 

910  GOSUB  1460 

920  CALL  "DIFS", A, A 

930  S1=S1»10*3/<H*9.8044) 

940  GOSUB  1400 
950  GOSUB  1460 
960  GO  TO  1230 
970  GOSUB  1340 
980  S1=S 
990  GOSUB  1460 
1000  PRINT  SD: 

1010  CALL  "INT",A,A 

1020  S1>S1«H 

1030  GOSUB  1280 

1040  GOSUB  1460 

1050  CALL  "DIF3”,A,A 

1060  CALL  "DIF3",A,A 

1070  S1=»S1»10^3/<H*2»9.8044) 

1080  GOSUB  1400 
1090  GOSUB  1460 
1100  GO  TO  1230 
1110  GOSUB  1400 
1120  S1=«S 
1130  GOSUB  1460 
1140  PRINT  SD:  '•»" 

1150  CALL  '•INT",A,A 

1160  S1»SHH  . 0E-3»H*9.8044 

1170  GOSUB  1340 

1180  GOSUB  1460 

1190  CALL  "INT",A,A 

1200  S1=S1»H 

1210  GOSUB  1280 

1220  GOSUB  1460 

1230  UIEUPORT  0,U3, 0,100 

1240  UINDOU  0,W3, 0,100 

1250  HOUE  @0:0,0 

1260  PRINT  @D:"»  oriq.  tiwe  record"; 
1270  RETURN 

1280  VIEWPORT  5,W3-30,4,32 
1290  V*»"OISPLACEMENT" 

1300  G««"D  nex  :” 

1310  M*«"D  Min  : " 

1320  N9-"  tnnl" 

1330  RETURN 

1340  VIEWPORT  5 , W3-30 , 36 , 64 
1350  V**"VELOCITY" 

1360  G««"V  nex  :" 

1370  M9«"V  Min 
1380  N««"  lM/»ec]" 


1390  RETURN 

1400  VIEWPORT  5,W3-30,68,96 
1410  V$="ACCELERATION" 

1420  G««"A  nax  ; 

1430  H**“A  Min  : ” 

1440  N*="  Cgl" 

1450  RETURN 

1460  CALL  "MAX", A, VI ,14 
1470  CALL  "MIN", A, V2, 13 
1480  WINDOW  0,N0,V2,V1 
1490  IF  D»32  THEN  1550 
1500  HOVE  8Di0,0 
1510  FOR  lal  TO  NO 
1520  DRAW  8D;I,A(I) 

1530  NEXT  I 

1540  GO  TO  1560 

1550  CALL  "DISP",A 

1560  AXIS  8D:400,0 

1570  FOR  I»0  TO  NO  STEP  2000 

1580  MOVE  eD;I,0 

1590  PRINT  eD:"HJ",I*H; 

1600  NEXT  I 

1610  PRINT  @D! "HHHHHHHHHHHIMSl" 

1620  MOVE  eD:4200,vT 

1630  PRINT  8Di  USING  1 670 : “ J " ,G* , VI »S1 , N* ; 

1640  MOVE  eD; 4200, VI 
1650  PRINT  eOi "J"  j 

1660  PRINT  8Di  USING  1 670 : " J “ ,M* , V2*S1 , Ni ; 

1670  IMAGE  A,7A,6D. 1D,11A 
1680  MOVE  8Di0,0.5«<V2-^Vl) 

1690  PRINT  9Di"HH"; 

1700  FOR  I»1  TO  LEN<V6)/2 
1710  PRINT  SD! "K  "; 

1720  NEXT  I 

1730  FOR  I-l  TO  LEN<V*) 

1740  X*-SEG<V6,I,1 ) 

1750  PRINT  SDiXf; "HJ"; 

1760  NEXT  I 
1770  RETURN 

1780  REM  PROCESS 

1790  CHARSIZE  3 
1800  CALL  "TIME",B« 

1810  PRINT  9Di  "1" 

1820  PRINT  e32,26;2 
1830  IF  F1<1  THEN  1890 
1840  IF  FI >1000  THEN  1890 
1850  IF  F2<1  THEN  1900 
1860  IF  F2>1000  THEN  1900 
1870  IF  F2<*F1  THEN  1900 
1880  GO  TO  1930 
1890  IF  C0=1  THEN  1930 

1900  PRI  "ANALYSIS  FREQUENCY  IS  OUT  OF  RANGE  ENTER  NEW  VALUES  OF  F1,F2G" 
1910  INPUT  F1,F2 
1920  CO  TO  1800 

1930  DELETE  F5,X5,V5,G5,T5, Y5,Z5,C,K ,A 
1940  DIM  A(N0) 

1950  C3-INT< (MEM0RY/8-18)/8) 

1960  DIM  F5<C0) ,X5<C0) ,V5<C0) ,C5(C0) ,T5(C0) ,Y5(C0) ,Z5(C0) ,C(C0) 

1970  FIND  B3!D1 

1980  READ  S3;W«,R2,H,T1 ,M,S 
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1990  FIND  eSiDl 

2000  E1=INT(1000/(K1»F8*H)) 

2010  E2*INT(  1000/(K2»F8t*H)  ) 

2020  E3=INT(1000/(K3»F8»H) ) 

2030  F5(l )=F1 

2040  IF  C0=1  THEN  2080 

2050  FOR  1=2  TO  CO 

2060  F5(I)  =  <F2/F1 ) * < 1/<C0-1  )  )*F5< I-l  ) 

2070  NEXT  I 
2080  X5=0.01 
2090  V5=0.01 
2100  05=0.01 
2110  T5=0.01 
2120  Y5=0.01 
2130  25=0.01 
2140  11=1 
2150  12=0 

2160  IF  Fl<10  THEN  2200 
2170  12=1 

2180  IF  FKIOO  THEN  2200 

2190  12=2 

2200  GO  TO  2370 

2210  E=E1 

2220  GQSUB  2460 

2230  IF  I1>C0  THEN  2310 

2240  GO  TO  2370 

2250  E=E2 

2260  GQSUB  2460 

2270  IF  I1>C0  THEN  2310 

2280  GO  TO  2370 

2290  E=E3 

2300  GOSUB  2460 

2310  DELETE  A 

2320  CALL  ’'TIME",T$ 

2330  GOSUB  3320 
2340  PRINT  @32,26:0 

2350  PRINT  eD:"_MAX.  POSSIBLE  ANALYSIS  FREQUENCIES 
2360  RETURN 
2370  F3=F5<I1) 

2380  C1=I1 

2390  12=12+1 

2400  FOR  11=1  TO  CO 

2410  IF  F5<I1)>10*I2  THEN  2430 

2420  NEXT  II 

2430  F4=F5( Il-l  ) 

2440  GO  TO  12  OF  2210,2250,2290 
2450  RETURN 

2460  READ  @3 : W* , R2 , H , T1 ,M , S , A2<M ,2 ) , A 

2470  FIND  e3:Dl 

2480  IF  E<=1  THEN  2550 

2490  J=1 

2500  FOR  1=1  TO  NO  STEP  E 
2510  A(J)=A(I> 

2520  J=J+1 
2530  NEXT  I 
2540  GO  TO  2560 
2550  E=1 
2560  N5»N0/E 
2570  H2»H»E 

2580  GO  TO  R2  OF  2590,2630,2660 


2590  CALL  "DIF3",A,A 

2600  CALL  "DIFS", A, A 

2610  S2=S 

2620  GO  TO  2670 

2630  CALL  "DIF3",A,A 

2640  S2=S»H2 

2650  GO  TO  2670 

2660  S2=S»10''-3*H2*2»9.8044 

2670  C2=I1-1 

2680  P1=0 . 0062831853«F3xH2 
2690  FOR  I=C1  TO  C2 
2700  C(I)=2»C0S(P1 ) 

2710  IF  C1=C2  THEN  2730 
2720  P1=P1*<F4/F3)*(1/<C2-C1)) 

2730  NEXT  I 

2740  IF  Pl<0.902  THEN  2790 
2750  IF  143.566/H2>F2  THEN  2790 

2760  PRINT  "LGUARNING:  FREQUENCIES  ARE  TOO  HIGH  FOR  YOUR  TIME  INCREMENT" 
2770  PRINT  "  ERRORS  ARE  OUTSIDE  THE  RANGE  -7  TO  +3  PERCENT  FOR" 

2780  PRINT  "  FREQUENCIES  ABOVE  " 5 143 . 566/H2 ; "  HERTZ" 

2790  FOR  I=C1  TO  C2 
2800  F=F5(I) 

2810  IF  F<=1  THEN  2840 
2820  N=1000/(H2»5*LGT<F) ) 

2830  IF  N<N5  THEN  2850 

2840  N=N5 

2850  X1=0 

2860  X2=0 

2a.’0  R3=0 

2880  R4»0 

2890  IF  F=0  THEN  2920 
2900  N4=N-l-500/<F»H2) 

2910  IF  N4>3  THEN  2930 
2920  N4»3 

2930  FOR  J*1  TO  N-1 
2940  X3=C(I)*X2-X1+A( J) 

2950  IF  J<N4  THEN  3000 

2960  IF  ABS(X3)<ABS<R4)  THEN  3030 

2970  R4-X3 

2980  T4aJ 

2990  GO  TO  3030 

3000  IF  ABS(X3)<AEiS(R3)  THEN  3030 

3010  R3=X3 

3020  T3=J 

3030  X1=X2 

3040  X2>X3 

3050  NEXT  J 

3060  IF  ABS(R3) >ABS(R4)  THEN  3090 

3070  R3=R4 

3080  T3=T4 

3090  X=S2»R3 

3100  V=0 . 00628318«F«X 

3110  G=6.288318*F»V/9,8044 

3120  T=Tl+(T3-0 .5)»H2 

3130  Y=0 . 00628318*F»S2»R4 

3140  Z=1 .OE-3*(T1+(T4-0.5)*H2)*F+0.25 

3150  Z=Z-INT(Z) 

3160  IF  Y>0  THEN  3180 
3170  Z=Z-0.5 

3180  IF  Z<*0.5  THEN  3200 
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3190  Z*Z-1 
3200  Z=360»Z 
3210  FSID’^F 
3220  X5(I)aX 
3230  y5(I>=»V 
3240  G5<I)»G 
3250  T5(1)*T 
3260  Y5(I)=Y 
3270  Z5(I)=Z 

3280  PRINT  @0;  USING  3290 : I ,F,V,N; 

3290  IMAGE  20T , 3D, 4X , 4D . 1 D , 4X , 5D . 2D,3X,6D 
3300  NEXT  I 
3310  RETURN 

3320  REM  ******  SUBROUTINE  TIME  ■»«•»•»•»*•»»«■»*»*»»*■»»*■««««•****•)(««•»•*«•*»**» 
3330  DIM  B<2) 

3340  C»»B« 

3350  FOR  1=1  TO  2 
3360  X«=SEG(C«,11 ,8) 

3370  Y*=SEG(X*,4,2) 

3380  Z4*SEG<X«,7,2) 

3390  B<I)®VAL(X*)+<UAL(Y$)»60+VAL<Z*> )/3600 
3400  C$»Ti 
3410  NEXT  I 
3420  U=B(2)-B(1) 

3430  B3=(U-INT<U) )*60 
3440  B4=<B3-INT<B3))*60 

3450  PRINT  SDi  “.^CALCULATION  TIME  (h^  ;  win  •.  sec  )  : 

3460  PRINT  8D!  USING  3470 : INT <U INT( B3 INT ( B4 ) 

3470  IMAGE  2D, A, 2D, A, 2D 
3480  RETURN 

5000  REM  - 

5010  REM  - 

5020  REM  THIS  PROCEDURE  MAY  BE  USED  TO  CALCULATE  SHOCK  SPECTRA  OF 

5030  REM  TIME  RECORDS  STORED  ON  EXT.  TAPE  <FILE  1  TO  7)  AND  TO  STORE 

5040  REM  THE  ACCORDING  DATA  ON  INT.  TAPE  (FILE  »  TO  FILE  »+7>, 

5050  FOR  Dl=l  TO  7 
5060  GOSUB  52 
5070  N=Ol+0 
5080  FIND  N 

5090  MARK  1 ,7000 

5100  FIND  N 

5110  WRITE  C0,Wt,F5,X5,V5,G5,T5,Y5,Z5 
5120  NEXT  D1 
5130  RETURN 
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